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Abstract 

Primary  producers  play  a  critical  role  in  the  oceanic  food  chain  and  the  global  cycling  of 
carbon.  The  marine  diazotroph  Trichodesmium  is  a  major  contributor  to  both  primary 
production  and  nitrogen  fixation  in  the  tropical  and  subtropical  oceans.  These  regions  are 
often  characterized  by  low  phosphorus  (P)  concentrations,  and  P  deficiency  of 
Trichodesmium  could  limit  growth,  and  constrain  nitrogen  fixation  by  this  genus.  Three 
key  components  of  P  metabolism  in  cyanobacteria  include:  inorganic  phosphate  (Pi) 
uptake,  dissolved  organic  phosphorus  (DOP)  hydrolysis,  and  polyphosphate  (polyP) 
biosynthesis  and  catabolism.  These  three  pathways  lay  the  foundation  for  the  four 
chapters  that  make  up  the  body  of  this  thesis.  Chapter  two  examines  the  relative 
importance  of  Pz  and  DOP  in  Trichodesmium  in  the  Sargasso  Sea  by  measuring  uptake 
rates  and  uptake  kinetics  of  ^^Pz  and  DOP  (using  alpha  labeled  adenosine-5 ’-triphosphate 
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P-ATP  as  a  model  P  ester  compound).  The  in  situ  uptake  rates  suggest  that  the 
contribution  of  P  ester  to  total  P  uptake  can  be  greater  than  25%  for  Trichodesmium,  and 
therefore  P  esters  are  important  for  supplementing  Trichodesmium’ s  P  demand.  Chapter 
three  confirms  the  presence  of  polyP  in  the  Sargasso  Sea  populations.  This  is  the  first 
detection  of  polyP  in  an  oligotrophic  system  and  may  indicate  micro-scale  variability  in 
DIP  or  DOP  supply.  Chapter  four  examines  the  regulation  of  several  genes  involved  in  P 
scavenging  in  cultures  of  Trichodesmium  erythraeum  IMS  101:  two  copies  of  a  high 
affinity  phosphate  binding  protein  (pstS  and  sphX)  and  two  putative  alkaline 
phosphatases  (phoA  and  phoX).  Chapter  five  calibrates  the  expression  of  phoX  to 
decreases  in  nitrogen  fixation  rate  in  Trichodesmium  cultures  allowing  the  use  of  phoX 
expression  as  a  molecular  indicator  of  P  limitation  of  nitrogen  fixation.  Detection  of 
phoX  expression  in  the  North  Atlantic  and  South  Pacific  identified  the  extent  and 
magnitude  of  P  limitation  in  these  regions.  Taken  together  this  thesis  provides  new 
insights  into  how  Trichodesmium  survives  in  low  P  environments,  utilizes  DOP,  stores 
polyP,  and  the  extent  to  which  P  can  limit  nitrogen  fixation  by  this  genus. 
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Chapter  One 


Introduction 

Phytoplankton  play  a  key  role  in  the  biogeochemistry  of  the  ocean  by  fixing 
carbon  dioxide  (CO2)  and  influencing  the  flux  of  carbon  to  the  deep  sea.  Carbon  export 
in  the  ocean  impacts  the  atmospheric  CO2  concentration,  and  therefore  climate,  over 
geological  time  scales  (Longhurst,  1991).  Therefore,  a  fundamental  question  in 
biological  oceanography  is  what  controls  primary  production  in  the  ocean.  In  many 
regions,  particularly  the  ocean  gyres,  the  concentration  of  nutrients  such  as  nitrogen  (N), 
phosphorus  (P)  and  iron  (Fe)  are  very  low,  and  therefore  they  are  thought  to  be  limiting 
to  primary  production.  Some  genera  of  cyanobacteria,  such  as  Trichodesmium,  are 
capable  of  converting  N2  gas  into  biologically  available  nitrogen  compounds  (termed  N2 
fixation).  Thus,  they  contribute  to  carbon  cycling  directly  (as  primary  producers)  and 
indirectly  by  fixing  atmospheric  N2,  which  can  in  turn  promote  the  growth  of  other 
phytoplankton  species  (Coles  et  ah,  2004).  Trichodesmium  populations  contribute  an 
estimated  8-47%  of  the  total  primary  productivity  in  the  tropical  North  Atlantic, 
(Carpenter  et  ah,  2004).  The  annual  input  of  “new  nitrogen”  by  Trichodesmium  has  been 
estimated  to  be  as  much  as  1.6  x  10  mol  N  into  the  tropical  North  Atlantic,  which  can 
equal  or  exceed  vertical  fluxes  of  nitrate  to  the  surface  (Capone  et  ah,  1997).  Therefore, 
this  genus  has  a  significant  influence  on  the  biogeochemistry  of  its  environment. 

Since  Trichodesmium  can  fix  atmospheric  N2  they  are  not  N-limited,  and  therefore 
the  productivity  of  this  organism  is  limited  by  other  physical  and  chemical  properties. 
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Because  of  its  biochemical  significance  there  have  been  a  number  of  studies  seeking  to 
identify  what  limits  Trichodesmium  and  N2  fixation.  In  particular,  Fe  (Berman-Frank  et 
ah,  2001;  Webb  et  ah,  2001;  Kustka  et  ah,  2002;  Moore  et  ah,  2009),  P  (Wu  et  ah,  2000; 
Sanudo-Wilhelmy  et  ah,  2001;  Dyhrman  et  ah,  2002;  Mulholland  et  ah,  2002;  Moutin  et 
ah,  2005;  Hynes  et  ah,  2009)  or  a  combination  of  the  two  (Mills  et  ah,  2004)  have  been 
repeatedly  identified  as  the  limiting  nutrients  for  Trichodesmium  in  various  regions  of  the 
ocean.  In  general,  these  studies  suggest  that  P  is  limiting  in  the  North  Atlantic  while  Fe 
is  probably  limiting  in  the  South  Atlantic  and  the  Pacific  (Moore  et  ah,  2004).  Mainly 
due  to  the  fact  that  during  stratified  periods  in  the  North  Atlantic  inorganic  P  (Pi) 
concentrations  can  be  two  orders  of  magnitude  lower  than  the  North  Pacific  (Cavendar- 
Bares  et  ah,  2001).  Meanwhile,  the  North  Atlantic  receives  a  large  influx  of  iron-rich 
dust  from  the  North  Africa,  and  therefore  has  higher  Fe  concentrations  than  the  Pacific. 
The  focus  of  this  thesis  is  to  better  understand  P  metabolism  in  Trichodesmium,  and 
better  evaluate  where  and  when  P  limits  Trichodesmium  growth  and  N2  fixation  in  the 
field. 

Three  key  pathways  that  mediate  P  metabolism  in  cyanobacteria  are  (1)  inorganic 
phosphate  uptake,  (2)  dissolved  organic  phosphorus  (DOP)  hydrolysis,  and  (3) 
polyphosphate  (polyP)  metabolism.  Fi  is  generally  considered  the  preferred  form  of  P 
because  it  is  directly  available  for  growth  (Vershinina  and  Znamenskaia,  2002).  In  Fi 
deficient  conditions,  some  microbes  can  scavenge  P  from  DOP  bond  classes  such  as  P 
esters  and  phosphonates.  In  the  upper  water  column  of  oligotrophic  regimes  DOP  often 
comprises  a  significant  portion  of  the  total  P  (Orrett  and  Karl,  1987;  Karl,  1999; 
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Ammerman  et  al.,  2003),  and  several  recent  studies  have  suggested  that  DOP  may 
support  from  12-60%  of  primary  productivity  in  the  North  Atlantic  (Mather  et  al.,  2008; 
Lomas  et  al.,  2009).  However,  it  is  difficult  to  chemically  characterize  the  composition 
of  the  DOP  pool,  and  therefore,  assaying  the  bio-availability  of  DOP  to  a  given  microbial 
group  remains  an  ongoing  challenge.  Chapter  two  of  this  thesis  examines  the  relative 
importance  of  Pz  and  DOP  to  supporting  growth  in  Trichodesmium  by  measuring  Pz  and 
DOP  (using  alpha  labeled  adenosine-5 ’-triphosphate  as  a  model)  uptake  and  uptake 
kinetics  in  the  Sargasso  Sea. 

If  P  uptake  is  in  excess  of  cellular  needs,  Pz  can  be  stored  as  polyP  granules. 
PolyP  is  a  linear  chain  of  phosphate  residues,  and  could  represent  an  important  reservoir 
for  P  when  this  nutrient  is  not  readily  available  in  the  environment.  Additionally,  polyP 
has  recently  been  identified  as  an  important  component  of  the  P  cycle  in  coastal  systems 
(Diaz  et  al.,  2008).  PolyP  has  previously  been  observed  in  Trichodesmium  populations 
from  the  Caribbean  Sea  (Romans  et  al.,  1994),  but  it  has  never  been  examined  in  a  low  P 
environment  such  as  the  Western  North  Atlantic.  Determining  how  Trichodesmium 
allocates  P  under  P  limiting  conditions  and  if  polyP  formation  occurs  will  provide  new 
insights  into  how  Trichodesmium  is  adapted  to  low  P  environments.  Chapter  three 
examines  polyP  formation  and  degradation  in  Trichodesmium  from  culture  and  from  the 
Sargasso  Sea. 

When  microbes  are  P-limited  they  induce  a  set  genes  whose  products  are  involved 
in  P  scavenging  including  genes  involved  in  Pz  uptake  and  DOP  hydrolysis.  Identifying 
which  genes  are  responsive  to  P  limitation  in  Trichodesmium  will  allow  for  a  better 
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understanding  of  how  Trichodesmium  senses  and  responds  changes  in  P  physiology.  Pi 
uptake  in  many  microbes  occurs  via  the  P  specific  transport  (Pst),  which  includes  the 
high  affinity  phosphate  binding  protein,  encoded  by  sphX  or  pstS  (Vershinina  and 
Znamenskaia,  2002;  Suzuki  et  ah,  2004).  DOP  hydrolysis  can  be  mediated  by  many 
different  enzymes,  including  alkaline  phosphatases  encoded  by  the  phoX  and  phoA  genes 
(Roy  et  ah,  1982;  Moore  et  ah,  2005;  Wu  et  ah,  2007).  Chapter  four  examines  the 
presence  and  regulation  of  pstS,  sphX,  phoA  and  phoX  in  cultures  of  Trichodesmium 
erythraeum  IMS  101  under  different  P  physiologies  to  determine  which  gene  transcripts 
are  responsive  to  changes  in  P  supply  and  the  time  scales  of  these  responses. 

Finally,  chapter  five  measures  the  extent  of  P  limitation  of  Trichodesmium  in  the 
North  Atlantic  and  South  Pacific.  Various  techniques  have  been  utilized  to  assess  algal 
physiology  in  situ.  The  results  of  these  experiments  can  vary  greatly  depending  on  the 
time  scale  considered  and  the  process  being  measured.  Nutrient  bioassays  are  probably 
the  most  common  technique  employed  in  the  field  to  assess  nutrient  stress;  in  this  method 
a  nutrient  (such  as  N  or  P)  is  added  to  phytoplankton  incubations,  and  community 
responses  such  as  chlorophyll  or  N2  fixation  rates  are  measured.  This  technique  is 
problematic  for  Trichodesmium  because  cell  lysis  generally  occurs  quickly  after 
Trichodesmium  colonies  are  placed  in  incubation  bottles  (Bums  et  ah,  2006). 
Additionally,  contamination  of  nutrients  and  bacteria  is  a  constant  problem  with  bottle 
incubations.  Finally,  when  marine  microbes  are  sampled  and  incubated  in  a  bottle  there 
tends  to  be  major  shifts  in  the  community  structure  that  makes  extrapolating  bottle 
incubation  results  back  to  the  field  difficult. 
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The  cellular  ratio  of  C:N:P  may  be  used  to  assess  nutrient  physiology. 
Phytoplankton  tend  to  adhere  to  a  preferred  ratio  of  106C:  16N:  P  (Redfield  ratio)  for 
growth  (Redfield,  1934).  If  the  ratio  of  C:P  falls  above  the  Redfield  ratio,  then  the  cells 
are  considered  to  be  P  deficient  (Redfield,  1958;  Bertilsson  et  ah,  2003;  Heldal  et  ah, 
2003).  Measures  of  cellular  stoichiometry  can  be  extended  to  include  trace  metals,  and  to 
infer  trace  metal  limitation  (Morel  et  ah,  1991;  Berman-Frank  et  ah,  2001).  However, 
cellular  stoichiometry  has  been  found  to  be  variable  in  Trichodesmium  spp.  (Krauk  et  ah, 
2006;  White  et  ah,  2006),  indicating  that  cellular  stoichiometry  is  not  fixed  and  non- 
Redfield  ratios  may  not  be  diagnostic  of  nutrient  limitation.  Another  commonly  used 
assay  of  nutrient  physiology  is  measuring  the  activity  of  enzymes  that  are  diagnostic  of 
particular  physiological  states.  For  example,  bulk  alkaline  phosphatase  activity  (APA)  is 
often  used  as  a  proxy  for  P  limitation  (Ammerman  and  Glover,  2000;  Dyhrman  and 
Ruttenberg,  2006).  Although  valuable,  none  of  these  techniques  are  cell-specific. 
Different  organisms  within  the  same  community,  or  even  different  species  within  the 
same  genus,  are  experiencing  and  are  adapted  to  different  microenvironments;  they  will, 
therefore,  have  different  physiologies  (Lomas  et  al.,  2004).  Thus,  key  details  may  be  lost 
by  averaging  cellular  responses  over  the  community  as  a  whole.  Alkaline  phosphatase 
activity  measured  using  enzyme-labeled  fluorescence  (ELF)  is  cell- specific,  but  is 
qualitative,  and  therefore  it  can  be  difficult  to  assess  relative  degrees  of  P  limitation 
(Dyhrman  et  al.,  2002).  Molecular  techniques  promise  to  be  a  useful  and  innovative 
approach  to  probing  cellular  physiology  in  a  species-specific  manner:  allowing 
researchers  to  distinguish  between  species-specific  and  community  responses.  Chapter 
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five  uses  the  expression  of,  phoX,  which  was  shown  to  be  regulated  by  P  physiology  in 
chapter  four,  as  a  marker  of  P  limitation.  Further,  the  regulation  of  this  gene  is  calibrated 
to  decreases  in  N2  fixation  rate  to  identify  regions  in  the  ocean  where  N2  fixation  is 
reduced  due  to  P  limitation. 

This  work  represents  a  detailed  study  of  P  metabolism  in  Trichodesmium.  Our 
approach  included  studies  of  Trichodesmium  physiology  in  the  field  and  culture  using 
gene  expression  patterns,  tracer-radio  isotope  experiments,  and  studies  of  cellular  P 
allocation.  Thus,  this  work  provides  new  insights  into  how  Trichodesmium  survives  in 
low  P  environments,  utilizes  DOP  versus  Pi,  stores  polyP,  and  the  extent  to  which  P  can 
limit  nitrogen  fixation  by  this  genus. 
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Abstract 


The  dynamics  and  relative  importance  of  inorganic  phosphate  (Pi)  and  dissolved 
organic  phosphorus  (DOP)  uptake  were  examined  in  Trichodesmium,  and  the  microbial 
community.  Trichodesmium  DOP  and  Fi  uptake  rates  were  also  compared  to  literature 
values  from  several  other  taxonomically  important  groups  in  the  Sargasso  Sea.  Uptake 

TO 

rates  and  uptake  kinetics  of  Pi  and  DOP  (using  alpha  labeled  adenosine-5 ’-triphosphate 

"JO 

P-ATP  as  a  model  P  ester  compound)  were  assayed  during  two  cruises.  The  in  situ 
uptake  rates  suggest  that  the  contribution  of  P  ester  to  total  P  uptake  can  be  greater  than 
25%  for  Trichodesmium,  the  microbial  community,  and  key  phytoplankton  groups 
(Prochlorococcus,  Synechococcus,  picoeukaryotes,  and  nanoeukaryotes),  depending  on 
the  location.  Based  on  the  kinetics  of  Pz  uptake  and  P  ester  uptake,  Trichodesmium  is  a 
poor  competitor  for  Pz  but  a  much  better  competitor  for  P  esters  as  compared  to  the 
microbial  community.  Trichodesmium  growth  rates  calculated  from  the  Pz  and  P  ester 
uptake  parameters  suggest  that  only  a  small  fraction  of  the  P  ester  pool  needs  to  be  bio- 
available  to  allow  reasonable  growth  rates.  These  data  underscore  the  importance  of  P 
esters  in  supporting  production  by  Trichodesmium  in  the  Sargasso  Sea,  and  provide  the 
first  in  situ  measurements  of  P  ester  uptake  and  uptake  kinetics  in  this  important  N2 
fixing  genus. 


Introduction 

Oceanic  gyres  represent  the  earth’s  largest  ecosystem,  and  have  a  significant 
effect  on  global  biogeochemical  cycles  and  climate  (Karl  2002).  Phosphorus  (P)  has 
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increasingly  been  identified  as  a  limiting  nutrient  in  the  ocean  gyres,  particularly  in  the 
Sargasso  Sea  region  of  the  North  Atlantic  (Wu  et  al.  2000;  Maranon  et  al.  2003). 
Furthermore,  it  is  P  rather  than  nitrogen  (N),  which  is  thought  to  be  the  limiting  nutrient 
over  geological  time  scales  due  to  inputs  of  N  from  N2  fixation  (Longhurst  1991). 
Dissolved  inorganic  phosphate  (Pi)  concentrations  in  gyres,  such  as  the  Sargasso  Sea,  are 
extremely  low,  ranging  from  sub-nanomolar  to  10  nmol  L'*  (Cavendar-Bares  et  al.  2001; 
Lomas  et  al.  2009),  and  numerous  approaches  have  identified  evidence  of  P  deficiency  in 
phytoplankton  from  this  region  (Wu  et  al.  2000;  Ammerman  et  al.  2003;  Lomas  et  al. 
2004). 

The  major  forms  of  P  in  the  upper  water  column  are  Pi  and  dissolved  organic 
phosphorus  (DOP).  Pi  is  considered  the  preferable  P  source  for  microbes  because  it  is 
directly  available  for  growth  (Vershinina  and  Znamenskaia  2002),  but  in  the  upper  water 
column  of  oligotrophic  regimes  DOP  often  comprises  a  significant  portion  of  the  total  P 
(Karl  and  Bjorkman  2001).  It  is  difficult  to  chemically  characterize  the  composition  of 
the  DOP  pool,  and  therefore,  assaying  the  bio-availability  of  DOP  to  a  given  microbial 
group  remains  an  ongoing  challenge.  Within  the  high  molecular  weight  DOP  pool  that  is 
amenable  to  characterization  (approximately  20-40%  of  the  DOP),  roughly  25%  has  a 
phosphonate  bond  and  roughly  75%  has  an  ester  bond  (Clark  et  al.  1998).  As  such,  P 
esters  represent  a  larger  potential  P  source  than  phosphonates,  and  are  the  focus  of  this 
study. 

One  enzyme  known  to  hydrolyze  P  esters,  alkaline  phosphatase,  has  received  a 
great  deal  of  attention  because  of  its  importance  to  microbial  DOP  hydrolysis  in  the 
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ocean  gyres,  and  because  there  are  commercially  available  substrates  to  measure  its 
activity  (Ammerman  and  Glover  2000;  Dyhrman  et  al.  2002;  Lomas  et  al.  2004).  Using 
community  alkaline  phosphatase  activity  (APA)  as  a  proxy  for  DOP  hydrolysis,  a  recent 
study  by  Mather  et  al.  (2008)  suggested  that  DOP  could  fuel  up  to  30%  of  the  primary 
productivity  in  the  North  Atlantic.  However,  both  phytoplankton  and  heterotrophic 
bacteria  contribute  to  community  hydrolysis  rates  and  cell-specific  assays  of  APA  have 
shown  that  there  can  be  a  great  deal  of  taxonomic  heterogeneity  in  the  activity  of  this 
enzyme  (Dyhrman  et  al.  2002;  Lomas  et  al.  2004;  Webb  et  al.  2007).  Both  of  these 
factors  would  influence  production  estimates  based  on  community  hydrolysis  rates. 

While  studies  measuring  APA  have  provided  insight  into  the  importance  and 
scale  of  P  ester  hydrolysis  in  the  ocean,  APA  assays  are  problematic  in  that  they  typically 
involve  adding  excess  substrate,  and  therefore  measure  a  maximal  rate  of  hydrolysis,  not 
the  in  situ  rate.  Another  approach  for  measuring  nutrient  utilization  involves  short-term 
incubations  with  tracer  levels  of  radio-labeled  compounds  to  measure  the  uptake  rate  of 
these  isotopes  into  cells  (Fu  et  al.  2005a;  Moutin  et  al.  2005;  Zubkov  et  al.  2007).  This 
method  has  the  added  advantage  of  allowing  direct  comparisons  of  both  the  uptake  of  Pi 
and  DOP,  and  has  been  previously  employed  in  a  culture  study  of  Synechococcus  (Fu  et 
al.  2006). 

There  have  been  several  studies  that  have  focused  on  Fi  uptake  in  cultures  (Fu  et 
al.  2005a)  and  field  populations  (Moutin  et  al.  2005;  Sohm  and  Capone  2006;  Sohm  et  al. 
2008)  of  Trichodesmium,  a  colony  forming  N2  fixing  cyanobacteria,  which  can  dominate 
the  phytoplankton  community  in  oceanic  gyres  such  as  the  Sargasso  Sea  (Capone  et  al. 
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1997;  Carpenter  et  al.  2004).  This  genus  has  received  much  attention  because  as  a 
diazotroph  its  growth  is  not  limited  by  N,  but  may  be  limited  by  P.  Additionally, 
Trichodesmium  can  be  an  important  source  of  N  to  other  co-occurring  microbes  (Orcutt 
et  al.  2001;  Mulholland  et  al.  2006).  In  previous  studies  of  Fi  uptake  it  has  been 
suggested  that  Trichodesmium  might  meet  a  majority  of  its  P  demand  from  DOP  in  the 
form  of  P  esters,  and  that  Trichodesmium  is  a  poor  competitor  for  Pz  relative  to  the  rest  of 
the  microbial  community  (Mulholland  et  al.  2002;  Sohm  and  Capone  2006).  However,  in 
situ  rates  of  DOP  uptake  have  not  been  measured  in  Trichodesmium  field  populations. 

In  this  study,  we  compare  the  variability  and  relative  importance  of  Pz  and  DOP 
uptake  in  Trichodesmium,  the  community,  and  several  taxonomically  important  groups. 
To  accomplish  this,  uptake  rates  and  uptake  kinetics  of  ^^Pz  and  DOP  (using  alpha  labeled 
adenosine-5’ -triphosphate  P-ATP  as  a  model  P  ester  compound)  were  assayed  during 
two  cruises  in  the  Sargasso  Sea. 

Methods 

Sample  collection 

Samples  were  collected  on  two  cruises  to  the  Sargasso  Sea  in  June  2007  (X0705) 
and  May  2008  (X0804)  (Fig.  1).  Trichodesmium  colonies  were  collected  from  the  near 
surface  (roughly  within  the  top  20  m)  using  a  hand  held  130  pm  net.  Single  colonies 
were  picked  into  0.2  /zm-filtered  local  surface  seawater.  Colonies  were  then  transferred 
into  fresh  0.2  /zm-filtered  water  to  reduce  contamination  of  closely  associated  organisms, 
and  were  separated  by  morphotype  (either  ‘puff’  with  radial  trichomes  or  ‘raft’  with 
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parallel  trichomes)  for  further  analyses.  Samples  of  the  bulk  microbial  community,  and 
for  flow  cytometric  sorting,  were  collected  from  a  Niskin  bottle  on  a  rosette  from  a  depth 
of  5  m  and  processed  as  described  below  and  in  Casey  et  al.  (2009). 

Uptake  measurements 

Uptake  measurements  for  Trichodesmium  and  the  bulk  microbial  community 
were  modeled  after  a  previous  Trichodesmium  uptake  study  (Sohm  and  Capone 
2006).  In  brief,  ten  washed  Trichodesmium  colonies  were  transferred  into  10  mL  of 
surface  seawater.  In  each  incubation  only  one  colony  morphoytype  (puff  or  raft)  was 
used  depending  on  which  morphotype  was  most  abundant  at  that  station.  To  account  for 
the  possibility  of  free  Trichodesmium  trichomes  contributing  to  the  microbial  community 
uptake  calculations,  4L  of  surface  seawater  was  filtered  onto  5/^m  filter,  and  counted.  At 
most  there  were  6  trichomes  L'\  and  therefore  it  is  unlikely  that  free  trichomes  had  a 
large  affect  on  uptake  rates  that  were  assayed  in  10  mL.  Coincident  with  the 
Trichodesmium  colony  sampling,  10  mL  of  surface  seawater  was  sampled,  and  both 
Trichodesmium  and  the  surface  seawater  (referred  to  as  the  microbial  community)  were 
incubated  with  either  37-93  kBq  of  H3^^P04  or  9.2-37  kBq  of  alpha  labeled  ^^P  -ATP  for 
20-40  or  40-70,  minutes  respectively.  Nucleotides  like  ATP,  were  the  only  commercially 
available  radio-labeled  P  esters  that  could  be  obtained  at  the  time  of  the  cruise  and  as 
such  P-ATP  was  used  as  the  model  for  P  ester  uptake.  ATP  has  previously  been  used  to 
examine  P  ester  uptake  in  cultures  of  Synechococcus  (Fu  et  al.  2006),  and  in  field 
populations  of  phytoplankton  (Casey  et  al.  2009),  also  allowing  direct  comparisons  to 
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these  studies.  The  H3  PO4  and  P  -ATP  were  sampled  in  sets,  and  the  differences  in 
incubation  time  reflect  the  time  needed  to  sample  one  set.  Data  are  referred  to  herein  as 

00  00 

Pi  or  P  ester  uptake  for  clarity.  Replicate  samples  for  both  Trichodesmium  and  the 
microbial  community  were  killed  with  1  mL  of  formalin  (1%  final  concentration)  to 
account  for  abiotic  adsorption  of  P.  The  activity  of  the  killed  control  (typically  about 
10%  of  the  total)  was  subtracted  from  the  total  measured  activity.  To  further  explore  the 
importance  of  adsorbed  P,  ^^Pi  uptake  was  assayed  as  above  on  an  oxalate  washed 
Trichodesmium  sample  and  a  non-washed  Trichodesmium  sample  following  the  protocol 
of  Sanudo-Wilhelmy  et  al.  (2004),  on  a  subsequent  cruise  in  a  similar  region  of  the 
Sargasso  Sea.  In  all  cases,  incubations  were  ended  by  filtration  onto  0.2  jim 
polycarbonate  filters,  and  then  5  mL  of  scintillation  cocktail  (AquaSol,  Perkin  Elmer) 
was  added  to  each  sample,  they  were  shaken,  and  counted  on  a  Tri-Carb  2900  liquid 
scintillation  counter  (Packard).  Assays  of  uptake  for  flow-sorted  populations  were 
performed  by  J.  Casey  and  co-authors  (Casey  et  al.  2009). 

Trichodesmium  uptake  rates  were  obtained  by  subtracting  the  microbial 
community  rates  from  the  Trichodesmium  and  microbial  community  combined  uptake 
values.  Samples  were  taken  in  triplicate,  but  on  occasion  Trichodesmium  biomass  was 
lost  to  the  walls  of  the  incubation  bottle  or  filter  tower,  and  one  of  the  replicates  had  to  be 
dropped.  Reported  uptake  rates  are  the  average  of  duplicate  or  triplicate  samples. 
Uptake  rates  were  calculated  using  the  following  formula:  Uptake  rate=  (Rs-Rk)  x  P  x  (A 
X  T  X  B)'^  where  Rs  is  the  radioactivity  of  the  sample,  Rk  is  the  radioactivity  of  the  killed 
control,  A  is  the  total  radioactivity  added  to  incubation,  T  is  the  time  of  incubation,  B  is 
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the  chlorophyll  a  concentration  of  a  replicate  sample  or  Trichodesmium  colony  number 
(used  for  growth  rate  calculations),  and  P  is  the  ambient  concentration  of  the  compound 

TO 

of  interest  in  the  sample.  For  P/  uptake  rates,  P  was  the  measured  soluble  reactive 
phosphorus  (SRP)  concentration.  The  ambient  ATP  concentration  was  not  measured  in 
this  study,  but  previous  bioassay  measurements  of  ATP  in  the  North  Atlantic  found 
concentrations  around  1  nmol  L''  (Zubkov  et  al.  2007).  Therefore,  an  ambient  ATP 
concentration  of  1  nmol  L'  was  used  to  estimate  a  minimum  P  ester  uptake,  as  has  been 
used  in  previous  work  (Casey  et  al.  2009).  However,  this  might  be  a  slight  overestimate 
of  this  minimum  value,  because  North  Pacific  ATP  concentrations  were  found  to  be  0.3- 
0.5  nmol  L'^  in  the  top  100  m  of  the  water  column  using  a  chemical  precipitation 
approach  (Karl  and  Bjorkman  2001).  A  maximal  P  ester  uptake  rate  was  calculated 
assuming  that  ATP  is  a  model  for  all  P  ester.  Since  75%  of  high  molecular  weight  DOP 
has  been  shown  to  be  P  ester  (Clark  et  al.  1998),  the  ^^P  ester  uptake  rates  are  reported  as 
a  range  using  both  1  nmol  L'^  and  75%  of  the  DOP  concentration  for  the  constant  P  in  the 
uptake  equation,  unless  otherwise  noted.  To  compare  uptake  rates,  all  values  were 
normalized  to  chlorophyll  a  after  Sohm  and  Capone  (2006).  However,  we  acknowledge 
that  there  was  heterotrophic  uptake  in  both  the  Trichodesmium  and  microbial  community 
samples,  as  has  been  noted  in  previous  studies  (Sohm  and  Capone  2006;  Sohm  et  al. 
2008). 

Uptake  Kinetics 
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For  each  kinetic  curve,  0-1.0  /^mol  L'^  cold  Pi  or  cold  ATP  was  added  in 
addition  to  the  Pi  or  P-ATP  spike  to  duplicate  incubations  and  assayed  as  above.  Data 
was  fit  (Prism,  Graphpad)  using  the  Michaelis-Menton  equation:  V=Vmax  x  S  (Km-i-S)'\ 
where  Vmax  is  the  maximal  uptake  velocity,  S  is  the  substrate  concentration  (either  the 
measured  Pi  concentration  or  the  estimated  value  of  1  nmol  L'^  for  P  ester),  and  Km  is  the 
nutrient  concentration  where  the  uptake  is  half  of  Vmax- 

Chlorophyll  a 

Ten  puff  or  raft  colonies  were  washed  and  filtered  onto  GF/F  filters  for 
Trichodesmium  analyses,  or  1-2  L  of  seawater  were  filtered  onto  a  GF/F  for  microbial 
community  analyses,  and  all  samples  were  frozen  at  -20°C.  Filters  were  extracted 
overnight  in  90%  acetone  at  -20°C.  Chlorophyll  a  (Chi  a)  fluorescence  was  measured  on 
an  Aquafluor  handheld  fluorometer  (Turner  Designs)  or  TD-700  fluorometer  (Turner 
Designs)  and  analyzed  according  to  Arar  and  Collins  (1997).  In  all  cases,  colonies  were 
chosen  to  best  mimic  the  size  of  colonies  used  for  the  other  assays  of  Trichodesmium,  and 
the  Chi  a  colony'^  is  reported  for  recalculating  the  normalized  data  as  necessary. 

P  Chemistry 

Soluble  reactive  phosphorus  (SRP)  was  measured  according  to  the  magnesium- 
induced  co-precipitation  method  (Karl  and  Tien  1992),  with  a  detection  limit  of  0.5  nmol 
L'\  SRP  concentrations  measured  here  are  referred  to  as  the  Pi  concentration  for 
consistency  with  the  reporting  of  the  uptake  data.  The  total  dissolved  P  (TDP)  was 


27 


processed  according  to  the  persulfate  oxidation  method  (Valderrama  1981).  DOP 
concentrations  were  calculated  as  the  difference  between  TDP  and  Pi. 

Alkaline  Phosphatase  Activity  (APA) 

For  Trichodesmium  APA  assays,  3-5  colonies  were  washed  and  placed  on  a  5  pm 
polycarbonate  filter.  For  microbial  community  APA,  300  mL  of  seawater  was  collected 
from  5  m,  filtered  onto  a  0.2  pm  filter  and  in  both  cases  the  filters  were  stored  frozen  at  - 
20‘’C.  Samples  were  processed  as  described  elsewhere  (Dyhrman  and  Ruttenberg  2006). 
Briefly,  10  /^mol  L'^  6,8-difluoro-4-methylumbelliferyl  phosphate  (di-MUF-P) 
(Invitrogen)  was  added  to  each  sample  in  a  petri  dish  with  artificial  seawater  containing 
no  added  P,  and  fluorescence  was  measured  on  a  Fluostar  Optima  plate  reader  (BMG 
Labtech)  every  5-20  minutes  for  5  time  points,  within  the  linear  range  of  the  assay. 
Previous  kinetics  experiments  found  the  10  /^mol  L’'  substrate  concentration  to  be 
saturating  (data  not  shown).  Standard  curves  were  generated  for  each  assay  using  6,8- 
difluoro-7-hydroxy-4-methylcoumarin  (di-MUF)  as  a  standard. 

Statistics 

Two-tail  paired  Mests  were  preformed  to  compare  data  sets.  Two  data  sets  were 
determined  to  be  significantly  different  from  each  other  if  they  had  a  p-value  less  than 
0.05. 
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Results 


The  Pi  concentrations  were  low  for  both  transects,  averaging  3.0  +3.1  nmol  L'^ 
(mean  +standard  deviation)  (Table  1).  The  DOP  concentration  was  more  than  15  times 
higher  than  the  Pi  concentration  with  an  average  of  48.0  +28.1  nmol  L'^  (p<0.001)  (Table 
1).  The  Chi  a  content  per  colony  was  similar  between  stations  with  an  average  of  28  +13 
ng  Chi  a  colony' '  as  was  the  microbial  community  Chi  a  concentration,  which  averaged 
31  +7.7  ng  Chi  a  L'^  (Table  1).  The  microbial  community  APA  averaged  60.7  +33.2 
nmol  P  h'^  jug  Chi  while  Trichodesmium  APA  had  an  average  of  19.5  +10.5  nmol  P  h' 
‘  //g  Chi  (Table  1). 

33  1 

The  average  Trichodesmium  Pz  uptake  for  both  cruises  was  0.98  +0.80  nmol  h' 
ng  Chi  d^  (Table  1).  Although  not  statistically  different,  the  average  microbial 
community  Pz  uptake  was  higher  than  that  of  the  Trichodesmium,  and  at  each  station 
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Pz  uptake  was  higher  for  the  microbial  community  than  for  Trichodesmium  (Table  1). 
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Furthermore,  the  average  proportion  of  Pz  uptake  by  Trichodesmium  (8%)  was 
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significantly  lower  (p<0.001)  than  the  average  proportion  of  Pz  uptake  by  the  microbial 
community  (92%)  relative  to  the  total  Pz  uptake  at  each  station.  Pz  can  adsorb  to 
Trichodesmium  colonies  and  other  constituents  of  the  microbial  community  (Sanudo- 
Wilhelmy  et  al.  2004).  To  account  for  this,  the  activity  assayed  in  killed  controls 
(approximately  10%  of  the  total),  was  subtracted  from  the  total  activity  to  calculate  the 
reported  uptake  rates.  Oxalate  washing  (Sanudo-Wilhelmy  et  al.  2004)  is  a  more 
stringent  way  to  remove  adsorbed  P.  A  comparison  of  Pz  uptake  on  an  oxalate  washed 
Trichodesmium  sample  vs.  a  non-washed  Trichodesmium  sample,  suggests  that  this  more 
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stringent  approach  results  in  26%  of  the  uptake  being  the  result  of  Pi  adsorbtion. 
Therefore,  the  values  presented  here  (calculated  with  10%)  might  be  a  roughly  16% 
overestimate  of  the  Trichodesmium  Pi  uptake  as  they  include  some  Pi  adsorbtion. 

TO 

Trichodesmium  rates  of  P  ester  uptake  were  variable,  with  an  average  minimum 
of  0.11+0.13  and  an  average  maximum  of  4.5  +7.1  nmol  P  h''/^g  Chi  using  1  nmol  L'^ 
or  75%  of  the  DOP  concentration  as  the  ambient  P  ester  concentration  respectively 
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(Table  1).  The  P  ester  uptake  rates  for  the  microbial  community  ranged  from  a 
minimum  of  1.4  +1.0  to  maximum  of  63  +65  nmol  P  h'  ng  Chi  calculated  with  ranges 
as  described  above  (Table  1).  Therefore,  for  both  Trichodesmium  and  the  microbial 
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community,  the  P  ester  uptake  ranged  from  lower  than  the  Pi  uptake  rate  to  higher, 
depending  on  what  the  ambient  bio-available  P  ester  concentration  is  assumed  to  be. 

To  compare  the  relative  importance  of  P  ester  vs.  Pz  uptake,  the  %P  ester  uptake 
of  the  total  P  uptake  (Pz  plus  P  ester)  for  Trichodesmium  and  the  microbial  community 
was  calculated  (Fig.  2).  For  Trichodesmium,  the  contribution  of  P  ester  to  total  P  uptake 
was  variable,  and  spanned  different  ranges  for  each  station  depending  on  the  ambient 
DOP  concentration  used  to  estimate  the  P  ester  concentration  (Fig.  2).  At  3  of  the  6 
stations,  P  ester  accounted  for  at  least  25%  of  the  P  uptake  in  Trichodesmium  and  the 
average  minimum  for  all  stations  was  29%  (Fig.  2).  P  ester  typically  contributed  a 
smaller  portion  to  total  P  uptake  for  the  microbial  community  (Fig.  2).  In  this  case,  P 
ester  accounted  for  at  least  25%  of  the  total  P  uptake  at  only  1  out  of  6  stations  for  the 
microbial  community,  and  the  average  minimum  contribution  was  10%  (Fig.  2). 
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Using  the  X0705  cruise  taxon  specific  uptake  rates  measured  by  Casey  et  al. 
(2009)  the  contribution  of  P  ester  to  total  P  uptake  was  calculated  for  picocyanobacteria, 
picoeukaryotes  and  nanoeukaryotes.  The  contribution  of  P  ester  to  total  P  uptake  in 
Prochlorococcus,  and  Synechococcus  was  variable  depending  on  the  station;  for  example 
at  several  stations,  P  ester  contributed  a  minimum  of  25%  to  the  total  P  uptake,  but  not  at 
all  stations  (Fig.  3).  In  contrast,  P  ester  contributed  a  minimum  of  25%  of  the  total  P 
uptake  for  every  station  assayed  for  the  picoeukaryotes  and  the  nanoeukaryotes.  At 
X0705  Sta.  16,  the  only  station  where  both  Trichodesmium  and  the  flow-sorted  data  are 
available,  the  contribution  of  P  ester  to  total  P  uptake  was  25%  or  more  for  each  group 
(Figs.  2,  3). 

'j'i 

Assays  for  Trichodesmium  and  the  microbial  community  Pi  uptake  kinetics  were 
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performed  at  four  stations;  one  station  on  X0705  and  three  stations  on  X0804.  P  ester 
uptake  kinetics  were  assayed  at  three  stations  on  cruise  X0804.  For  the  Pi  uptake 
kinetics,  Trichodesmium  had  an  average  Vmax  of  210  +120  nmol  P  h'^g  Chi  and  an 
average  Km  of  0.68  +0.33  //mol  L'*  (Table  2).  The  average  Km  value  is  similar  to  those 
obtained  in  previous  studies  with  non-oxalate  washed  samples  which  include  some 
adsorbed  P  (Table  3),  however  the  Vmax  calculated  herein  is  higher  than  those  previous 
studies  (Table  3).  The  Vmax  for  ^^P  ester  uptake  had  an  average  of  54  +34  nmol  P  h''  //g 
Chi  a'^  and  the  Km  averaged  0.29  +0.18  //mol  L'*  (Table  2).  For  the  microbial 
community  the  Km  and  Vmax  of  ^^Pi  uptake  averaged  0.057  +0.081  //mol  L'^  and  17  +9.1 
nmol  P  h'*//g  Chi  respectively  (Table  2).  The  Km  and  Vmax  of  ^^P  ester  uptake  for  the 
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microbial  community  averaged  0.51  +0.42  /rmol  L'^  and  11  +6.1  nmol  h'^  //g  Chi 
respectively,  in  the  microbial  community  samples  (Table  2). 


One  measure  of  maximal  P  ester  hydrolysis  and  uptake  is  the  V^ax  of  ester 
uptake  as  described  above.  An  alternative  measure  of  maximal  P  ester  hydrolysis  is 
APA.  The  P  ester  Vmax  for  Trichodesmium  was  similar  to  (within  the  same  order  of 
magnitude),  but  always  slightly  higher  than  the  APA  (Table  2).  Likewise,  ^^P  ester  Vmax 
for  the  microbial  community  was  also  similar  (within  the  same  order  of  magnitude)  to  the 
APA.  There  were  very  few  stations  where  both  the  APA  and  the  Vmax  of  P  ester  uptake 
were  measured  (Table  2),  therefore  none  of  the  kinetic  parameters  were  statistically 
significantly  different  from  each  other  due  to  the  limited  dataset. 


Discussion 

The  Sargasso  Sea  has  previously  been  identified  as  a  region  where  Pi  is  limiting 
to  phytoplankton  growth  (Wu  et  al.  2000;  Mather  et  al.  2008).  In  this  region,  being  able 
to  utilize  the  DOP  pool  as  a  P  source  may  be  an  important  growth  advantage  for 
microbes.  However,  determining  the  bio-availability  of  DOP  is  a  difficult  question  to 
address,  as  the  bio-availability  of  DOP  varies  among,  and  sometimes  within,  a  microbial 
taxa  (Lomas  et  al.  2004;  Moore  et  al.  2005).  Therefore,  determining  the  importance  of 
DOP  relative  to  Pi  uptake  remains  a  challenge,  but  is  central  to  understanding  how  P 
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controls  growth  in  this  region  of  the  ocean.  By  combining  measurements  of  Pi  and  P 
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ester  uptake  and  uptake  kinetics  and  we  can  gain  new  insights  into  the  importance  of 
DOP  to  Trichodesmium,  and  the  microbial  community,  in  the  low  P  Sargasso  Sea. 


Uptake  Rates 

The  Chi  a  normalized  Pi  uptake  rate  for  Trichodesmium  and  the  microbial 
community  were  variable  across  both  transects,  with  no  clear  trend  with  ambient  Fi  or 
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DOP  concentration.  However,  at  each  individual  station  the  Pz  uptake  rate  for 
Trichodesmium  was  at  least  an  order  of  magnitude  lower  than  that  of  the  microbial 
community,  and  the  proportion  of  total  Pz  uptake  for  Trichodesmium  was  significantly 
lower  than  the  community.  Therefore,  Trichodesmium  is  not  taking  up  as  much  Pz  per 
unit  Chi  a  as  the  microbial  community,  and  this  suggests  that  Pz  may  not  be  as  important 
a  P  source  for  Trichodesmium  as  it  is  for  the  microbial  community.  These  data  are 
consistent  with  previous  comparisons  of  Trichodesmium  and  the  microbial  community  in 
this  region,  which  suggest  that  Trichodesmium  is  a  poor  competitor  for  Pz  (Sohm  and 
Capone  2006).  It  is  important  to  note  that  heterotrophic  bacteria,  which  do  not  contain 
Chi  a,  are  present  in  the  microbial  community  and  associated  with  Trichodesmium,  and 
therefore  likely  contribute  to  measured  uptake  rates.  As  a  result,  the  assay  and 
normalization  approach  used  here  and  elsewhere  (Sohm  and  Capone  2006;  Sohm  et  al. 
2008)  may  overestimate  the  uptake  rates  of  Pz  or  P  ester  in  the  microbial  community,  and 
Trichodesmium.  Trichodesmium,  in  particular,  is  known  to  have  tightly  associated 
epibionts.  It  is  not  clear,  however,  whether  the  P  taken  up  or  hydrolyzed  by  epibionts 
would  eventually  contribute  to  meeting  the  P  demand  of  Trichodesmium. 


33 


Examining  the  percent  of  P  ester  uptake  contribution  to  the  total  P  uptake,  allows 
the  relative  importance  of  P  esters  as  a  P  source  to  be  assessed  for  Trichodesmium,  the 
microbial  community,  and  different  microbial  taxa.  This  calculation  also  cancels  out  the 
Chi  a  normalization  and  thus  eliminates  this  confounding  factor  in  comparing  the 
microbial  community  to  Trichodesmium.  Given  the  uncertainty  regarding  the  bio- 
available  P  ester  concentration,  this  proportion  was  calculated  to  determine  the  potential 
minimum  and  maximum  contribution  of  P  ester  to  total  P  uptake.  Although  it  is  known 
that  some  of  the  microbial  populations  in  this  region  have  the  genetic  machinery  to  take 
up  phosphonate  DOP  (Dyhrman  et  al.  2006;  Ilikchyan  et  al.  2009),  phosphonate  uptake 
should  be  smaller  than  combined  Pz  and  P  ester  uptake  as  these  chemical  pools  are  more 
abundant  and  considered  more  bio-available.  At  half  of  the  stations,  P  ester  contributed  a 
minimum  of  25%  of  the  total  P  uptake  for  Trichodesmium.  In  contrast,  there  was  only 
one  station  where  the  P  ester  contribution  was  greater  than  25%  of  the  total  P  uptake  for 
the  microbial  community.  These  data  underestimate  the  importance  of  P  ester  to 
Trichodesmium,  because  accounting  for  Pz  that  is  adsorbed  to  Trichodesmium  would 
reduce  Pz  uptake  rates  by  roughly  16%  but  should  not  affect  P  ester  uptake  rates  (Fu  et  al. 
2005b;  Fu  et  al.  2006).  Taken  together  these  results  suggest  that  P  esters  are  a  major  P 
source  for  Trichodesmium,  but  probably  less  important  to  the  microbial  community. 
Although,  for  both  Trichodesmium  and  the  microbial  community,  the  specific 
contribution  of  P  ester  to  total  P  uptake  is  highly  dependent  on  the  actual  concentration  of 
P  ester  that  is  bio-available. 
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Using  flow  cytometry  coupled  with  isotope  incubations  Casey  et  al.  (2009) 
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measured  Pi  and  P-ATP  (as  a  proxy  for  P  ester)  uptake  in  several  other  taxa  including 
Prochlorococcus,  Synechococcus,  picoeukaryotes  and  nanoeukaryotes.  Strikingly,  even 
at  the  average  lower  limit  of  P  ester  contribution,  P  ester  is  the  major  (greater  than  an 
average  of  50%  at  all  stations)  P  source  for  all  four  of  these  groups.  These  data  suggest, 
that  like  Trichodesmium,  these  critical  phytoplankton  groups  heavily  rely  on  P  esters  to 
support  P  demand  in  this  system.  This  is  consistent  with  culture  work  examining  ATP 
uptake  in  Synechococcus  that  suggests  that,  at  least  for  some  strains,  P  ester  is  important 
to  meeting  P  demand  (Fu  et  al.  2006).  At  Sta.  16  on  X0705,  there  were  simultaneous 
measurements  of  Trichodesmium,  microbial  community,  and  taxon  specific  Pi  and  P 
ester  uptake  allowing  for  direct  comparisons.  For  this  station,  P  esters  were  particularly 
important  for  all  individual  taxa,  even  at  the  lower  limit  of  P  ester  contribution,  although, 
the  microbial  community  does  not  reflect  this  trend.  This  observation  emphasizes  the 
fact  that  other  groups  such  as  heterotrophic  bacteria  are  probably  contributing  to  the 
microbial  community  uptake  measurements  and  are  good  competitors  for  Pi. 
Interestingly,  Sta.  16  is  one  of  the  stations  with  the  lowest  Fi  concentrations  of  <0.5  nmol 
L'\  This  suggests  that,  as  one  might  predict,  with  lower  Pi  concentrations  P  esters 
become  more  important  to  P  uptake. 


Kinetics 

The  average  Trichodesmium  Km  for  Pz  and  P  ester  uptake  were  similar.  Although 
both  affinity  constants  are  higher  than  the  ambient  P  ester  and  the  Pz  concentrations  of  the 
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upper  water  column  in  this  system,  the  P/  Km  is  consistent  with  the  Km  measured  in 
previous  culture  and  field  studies  which  ranges  between  179-710  nmol  L'^  for  assays  that 
did  not  remove  adsorbed  P  (Table  3).  These  may  be  overestimates  because  the  Km  of  Fi 
uptake  is  reduced  when  surface  adsorbed  P  is  removed  (Fu  et  al.  2005b).  Nevertheless, 
the  Km  values  detected  in  the  field,  across  a  range  of  Fi  concentrations,  are  consistent 
with  previous  work,  which  found  that  Trichodesmium  does  not  alter  its  Km  in  response  to 
P  physiology,  or  ambient  P  supply  (Fu  et  al.  2005a).  Given  that  the  concentration  of 
DOP  (and  thus  likely  bio-available  P  ester  concentration)  is  higher  than  that  of  Fi  in  the 
Sargasso  Sea,  Trichodesmium  is  potentially  better  at  taking  up  P  ester  than  Fi. 

Unlike  the  Km,  Trichodesmium  increases  its  ^^Pi  Vmax  in  response  to  increasing  P 
limitation  (Fu  et  al.  2005a).  In  this  study,  the  average  Vmax  of  Pi  uptake  for 
Trichodesmium  was  nearly  an  order  of  magnitude  higher  than  those  calculated  in  other 
studies  of  Trichodesmium  from  the  North  Atlantic  (Sohm  and  Capone  2006;  Sohm  et  al. 
2008)  which  also  did  not  remove  adsorbed  P.  There  are  two  additional  studies  that 
measured  Vspedfic  of  Trichodesmium  Pz  uptake:  Fu  et  al.  (2005a)  in  culture  and  Moutin 
et  al.  (2005)  in  the  South  Pacific.  Approximating  Vmax  from  those  studies  (Table  3)  also 
results  in  lower  Vmax  values  for  Trichodesmium  from  the  South  Pacific  (Moutin  et  al. 
2005),  and  from  P  replete  cultures  (Fu  et  al.  2005a)  relative  to  the  Vmax  for  ^^Pz  uptake 
calculated  here.  However,  the  average  Pz  concentration  from  the  stations  where  uptake 
kinetics  were  measured  in  this  study  was  much  lower  than  that  of  any  previous  field 
study  (4.6  nmol  L'*  vs.  45  nmol  L'^  for  the  second  lowest  Pz  concentration;  Table  3). 
Therefore,  the  higher  Vmax  values  found  in  this  study  may  reflect  Trichodesmium 
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populations  that  are  more  P- limited  than  those  in  previous  studies.  Fu  et  al.  (2005b)  found 
that  when  surface  adsorbed  P  was  removed  this  reduced  the  V^ax  for  Pz  uptake  by  half  for 
P-limited  Trichodesmium  cultures,  yet  a  50%  reduction  in  the  average  ^^Pz  Vmax 
calculated  herein  would  still  be  higher  than  other  field  studies  (Table  3).  This 
emphasizes  the  fact  that  the  Trichodesmium  in  this  study  are  in  a  particularly  low  Pz 
system,  and  may  be  P-limited. 

Comparisons  of  kinetic  parameters  between  Trichodesmium  and  the  microbial 
community  are  difficult  because  of  the  potential  for  differential  contributions  of 
individual  taxa  and  heterotrophic  groups.  While  the  Km  of  Pz  uptake  for  the  microbial 
community  is  highly  variable,  the  average  is  an  order  of  magnitude  lower  than  that  of 
Trichodesmium.  As  was  discussed  previously,  this  suggests  that  Trichodesmium  is  a  poor 
competitor  for  Pz  as  compared  to  the  microbial  community.  The  Km  of  P  ester  uptake  for 
the  microbial  community  is  similarly  variable,  but  the  average  is  higher  than  the  Km  of  P 
ester  for  Trichodesmium.  This  suggests  that  Trichodesmium  is  a  better  competitor 
compared  the  microbial  community  for  P  ester  than  for  Pz.  This  conclusion  is  consistent 
with  the  observations  with  Trichodesmium  made  above,  and  those  of  others  made  with 
assays  of  ^^Pz  alone  (Sohm  and  Capone  2006;  Sohm  et  al.  2008).  These  conclusions  are 
also  consistent  with  the  fact  that  P  ester  represented  25%  or  more  of  the  total  P  uptake  for 
Trichodesmium  at  more  stations  than  the  microbial  community. 

In  this  study,  the  maximal  rate  of  P  ester  utilization  was  examined  in  two  ways; 

TO 

through  P  ester  uptake  kinetics  and  through  APA  assays.  APA  measurements  here 
represent  the  maximal  rate  of  P  ester  hydrolysis  using  a  model  P  ester.  Whereas,  the  Vmax 
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of  P  ester  uptake  measures  the  maximal  rate  of  P  ester  hydrolysis  plus  incorporation 
into  the  cell,  assuming  that  the  P  ester  is  not  directly  taken  up  by  the  cell.  Each  of  these 
assays  uses  very  different  model  P  esters.  Given  the  inherent  differences  between  these 
approaches,  the  resulting  values  were  strikingly  similar  for  both  Trichodesmium  and  the 
microbial  community.  Together,  these  data  are  suggestive  of  a  population  that  can  take 
advantage  of  fluctuations  in  the  DOP  pool.  But  because  maximal  hydrolysis  rates  are 
always  much  higher  than  in  situ  rates,  these  data  underscore  the  fact  that  tracer  studies  are 
excellent  compliments  to  uptake  and  hydrolysis  rates  calculated  based  on  APA  (non¬ 
tracer)  alone. 


Growth  rate 

Determining  growth  rates,  or  doubling  times,  for  phytoplankton  provides  insights 
into  how  an  organism  is  adapted  to  its  environment,  and  its  nutrient  scavenging 
strategies.  Previous  studies  have  used  C  and  N  uptake  in  the  field  to  estimate  a 
Trichodesmium  doubling  time  (Carpenter  and  Romans  1991;  Orcutt  et  al.  2001).  These 
doubling  times  are  variable;  literature  values  range  from  0.7-2,500  days  (Orcutt  et  al. 
2001;  Moutin  et  al.  2005).  An  estimate  for  a  rapid  doubling  time  has  been  made  using  C 
uptake  of  roughly  2  days  for  Trichodesmium  in  the  Sargasso  Sea  (Orcutt  et  al.  2001). 
This  is  similar  to  doubling  times  for  Trichodesmium  in  culture  of  2-4  days  (Capone  et  al. 
1997;  Krauk  et  al.  2006).  Using  the  average  kinetic  parameters  from  our  study,  different 
P-specific  doubling  times  can  be  determined  at  a  range  of  Pz  and  P  ester  concentrations 
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by  following  the  equation  of  Moutin  et  al.  (2005)  (td)=  ln(2)(Km+S)(24  x  Vsp  x  S)'\ 
Where  td  is  the  doubling  time,  Vsp  is  the  V^ax  divided  by  the  P  quota  per  colony;  here  the 
P  quota  was  assumed  to  be  2.6  nmol  P  colony'^  (Sanudo-Wilhelmy  et  al.  2004),  and  S  is 
the  substrate  concentration  (Table  4).  It  is  important  to  note  that  the  P  quotas  are  from 
colonies  without  adsorbed  P  and  are  therefore  a  conservative  estimate  (Sanudo-Wilhelmy 
et  al.  2004).  Three  different  Pi  concentrations  were  used;  6  nmol  L''  (the  average  Pi 
concentration  for  the  stations  where  kinetic  parameters  were  measured),  3.3,  and  0.42 
nmol  L'^  (based  on  Zubkov  et  al.  (2007)  which  suggested  that  only  7-55%  of  the  Pi 
concentration  is  bio-available).  For  both  the  6  and  3.3  nmol  L'^  Pi  concentration,  the 
doubling  times  are  reasonable  for  the  Sargasso  Sea  and  roughly  equivalent  to  the 
Trichodesmium  doubling  times  calculated  using  C  (Orcutt  et  al.  2001).  For  the  P  ester 
calculations,  concentrations  of  1,  10,  and  36  nmol  L'*  were  used  based  on  the  possible 
range  of  bio-available  P  ester.  If  only  10  nmol  L"'  P  ester  is  bio-available,  then 
Trichodesmium  can  reach  the  same  doubling  time  as  on  6  nmol  L'^  of  Pi.  Therefore,  only 
a  small  portion  of  the  DOP  pool  needs  to  be  bio-available  P  ester  to  fully  support 
Trichodesmium  P  uptake,  based  on  these  data.  This  is  true  even  if  the  P  quota  is 
increased,  and  using  a  longer  doubling  time  in  this  analysis  only  reinforces  the  relative 
importance  of  P  ester  compared  to  Pi  in  this  system. 

Here  we  address  the  role  of  Pz  and  P  ester  in  supporting  the  growth  of 
Trichodesmium,  the  microbial  community,  and  draw  comparisons  with  different 
phytoplankton  groups.  The  in  situ  uptake  rates  suggest  that  the  contribution  of  P  ester  to 
total  P  uptake  can  be  greater  than  25%  for  Trichodesmium,  the  microbial  community,  and 
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key  phytoplankton  groups,  depending  on  the  station.  The  P  ester  contribution  to 
Trichodesmium  was  over  25%  at  half  of  the  stations  examined,  and  further  work  in  other 
systems  with  variable  P  chemistry  would  confirm  if  this  frequency  is  similar  in  higher  P 
systems,  or  to  what  extent  the  P  ester  uptake  rate  changes  with  ambient  DOP  over  a 
larger  dataset.  Several  studies  have  suggested  that  Trichodesmium  is  P-limited  in  the 
Sargasso  Sea  (Sanudo-Wilhelmy  et  al.  2001;  Dyhrman  et  al.  2002;  Mulholland  et  al. 
2002)  and  the  high  Fi  Vmax  values  found  here  support  that  conclusion.  A  comparison  of 
the  Km  values  for  both  P;  and  P  ester  uptake  suggests  that  Trichodesmium  is  a  better 
competitor  for  P  ester  than  for  Fi  relative  to  the  microbial  community.  This  is  consistent 
with  previous  work  (Sohm  and  Capone  2006),  which  focused  solely  on  Pz  uptake 
kinetics.  Future  work  focused  on  the  uptake  kinetics  of  other  individual  phytoplankton 
taxa  would  help  evaluate  the  consistency  of  this  relationship.  Additionally,  this  work 
focused  on  enzyme  activity  and  uptake  assays,  highlights  the  importance  of  pairing  tracer 
level  experiments  (which  measure  in  situ  uptake  rates)  to  experiments  where  saturating 
concentrations  of  substrate  are  added  (which  measure  Vmax  of  enzymatic  hydrolysis). 

-lO  O'! 

Trichodesmium  growth  rates  calculated  from  the  Pz  and  P  ester  uptake  parameters 
indicate  that  only  a  small  fraction  of  the  P  ester  pool  needs  to  be  bio-available  to  meet  P 
demand,  which  further  supports  the  importance  of  P  ester  to  this  group  in  the  Sargasso 
Sea.  Taken  together,  these  data  underscore  the  importance  of  P  esters  to  supporting 
production  by  Trichodesmium,  picocyanobacteria,  and  small  eukaryotes  in  the  Sargasso 
Sea. 
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Table  2.  Kinetic  parameters  of  P/  and  P  ester  uptake  assayed  for  Trichodesmium  and  the  bulk  microbial  community. 
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Table  3.  A  comparison  of  Trichodesmium  Vi  uptake  kinetics  and  P/  concentrations  in  this 
and  other  studies  of  Trichodesmium  from  cultures  and  field  populations. 


Source 

Location 

Pi 

(nmol  L ') 

Colony 

morphotype 

(nmol  L  ') 

33pj  V 

^  y  max 

(nmol  h  '  //g  Chi 

puffs  or 

This  study 

N.  Atlantic 

4.6 

rafts 

680 

210 

(Sohm  and  Capone  2006) 

N.  Atlantic 

65 

mixed 

710 

10.4 

(Sohm  et  al.  2008) 

N.  Atlantic 

45 

mixed 

180 

7.0 

(Moutin  et  al.  2005) 

S.  Pacific 

20-240 

mixed 

630 

3.1* 

P-replete 

NA 

340 

36* 

(Fu  et  al.  2005a) 

Cultures 

P-limited 

NA 

650 

80 

NA  -  not  applicable 

*  These  values  are  derived  assuming  a  P  quota  of  2.6  nmol  per  colony  from  Sanudo- 
Wilhelmy  et  al.  (2004)  and  an  average  Chi  a  content  per  colony  of  27  ng  (Table  1). 
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Table  4.  Trichodesmium  doubling  times  calculated  using  the  average  P/  and  P  ester 
kinetic  parameters  from  Table  2,  and  a  range  of  potentially  bio-available  P 
concentrations. 


P/  P  ester  Doubling  time 

(nmol  L'^) _ (nmol  L'^) _ (days) 

6  ^  Ta 

3.3  -  2.5 

0.42  -  20 

36  0.36 

10  1.4 

1  14 
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Figure  1.  Station  locations  for  two  cruises  to  the  Sargasso  Sea:  X0705  in  June  2007  and 
X0804  in  May  2008. 
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Figure  2.  The  estimated  range  of  P  ester  uptake  relative  to  total  P  uptake  for  (A) 
Trichodesmium  and  (B)  the  microbial  community.  The  black  bars  span  the  upper  and 
lower  limits  of  P  ester  uptake,  and  were  calculated  using  a  lower  limit  of  1  nmol  L'^  (the 
predicted  ATP  concentration)  and  an  upper  limit  of  75%  of  the  ambient  DOP 
concentration  (representing  the  potential  P  ester  concentration).  The  dashed  line 
indicates  the  average  lower  limit  of  P  ester  uptake  relative  to  total  P  uptake. 
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Figure  3.  The  estimated  range  of  the  proportion  of  P  ester  uptake  as  a  function  of  total  P 
uptake  for  Prochlorococcus,  Synechococcus,  picoeukaryotes,  and  nanoeukaryotes  during 
X0705  recalculated  from  Casey  et  al.  (2009).  The  black  bars  span  the  upper  and  lower 
limits  of  P  ester  uptake,  and  were  calculated  using  a  lower  limit  of  1  nmol  L'^  (the 
predicted  ATP  concentration)  and  an  upper  limit  of  75%  of  the  ambient  DOP 
concentration  (representing  the  potential  P  ester  concentration).  The  dashed  line 
indicates  the  average  lower  limit  of  P  ester  uptake  relative  to  total  P  uptake. 
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Abstract 


Polyphosphate  (polyP),  linear  ehains  of  three  or  more  phosphate  residues,  is  often 
considered  to  be  the  product  of  luxury  uptake  in  areas  of  excess  phosphorus  (P). 
However,  polyP  can  also  accumulate  in  P-depleted  cells  in  response  to  pulsed  P  re¬ 
supply.  As  such,  the  ability  to  synthesize  and  degrade  polyP  may  be  a  key  adaptation  to 
survive  in  environments,  like  the  oligotrophic  oceans,  with  a  variable  P  supply.  To  test 
the  hypothesis  that  polyP  is  present  in  phytoplankton  from  oligotrophic  systems,  the 
marine  diazotroph,  Trichodesmium  was  collected  from  the  low  P  surface  waters  of  the 
Sargasso  Sea  and  assayed  with  solid-state  ^^P  nuclear  magnetic  resonance  (NMR) 
spectroscopy.  In  this  system,  up  to  25%  of  Trichodesmium  cellular  phosphorus  was 
characterized  as  polyP,  despite  physiological  data  that  indicated  the  colonies  were  P 
deplete.  These  data  are  consistent  with  culture  studies,  which  underscore  that  all 
Trichodesmium  species  have  the  genetic  machinery  to  produce  polyP,  and  that  allocation 
of  P  to  polyP  increases  under  P  deplete  conditions  relative  to  P  replete  conditions.  Trends 
in  the  amount  of  Trichodesmium  polyP  along  the  cruise  transect,  show  that  allocation  of 
P  to  polyP  is  consistently  high,  and  that  the  ratio  of  polyP:C  varies  in  response  to  changes 
in  temperature  and  mixed  layer  depth.  We  hypothesize  that  Trichodesmium  is  taking 
advantage  of  fluxes  in  P  supply,  and  that  polyP  is  a  physiological  fingerprint  of  this 
variability.  Additionally,  if  polyP  formation  is  a  common  trait  in  phytoplankton,  these 
data  suggest  that  polyP  released  from  cells  is  an  additional  bio-available  component  of 
the  dissolved  organic  phosphorus  pool.  Taken  together,  this  study  highlights  the 
importance  of  polyP  to  P  cycling  and  cellular  P  allocation  even  in  oligotrophic  regions. 
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Introduction 


Inorganic  polyphosphate  (polyP)  is  a  polymer  of  phosphate  ranging  in  length 
from  3-1000’s  of  groups.  PolyP  synthesis  and  catabolism  are  typieally  controlled  by  a 
gene  {ppK)  eneoding  a  polyphosphate  kinase  that  reversibly  adds  phosphate  to  the  end  of 
the  polyP  ehain  (Tzeng  and  Kornberg  1998),  and  a  gene  {ppX)  encoding  an 
exopolyphosphatase  that  removes  the  terminal  phosphate  from  a  polyP  moleeule 
(Akiyama  et  al.  1993).  PolyP  has  been  found  in  all  major  groups  of  life  examined  to 
date,  but  its  funetion  is  varied,  and  in  many  regards  remains  unelear  (Kornberg  et  al. 

1999) .  Aceumulation  of  eellular  polyP  has  been  variously  attributed  to  a  stationary  phase 
adaptation,  an  energy  storage  eompound,  an  osmotie  regulator,  a  buffer  against  alkali 
conditions,  a  factor  in  DNA  competency  (as  part  of  a  DNA  channel),  and  in  intracellular 
phosphorus  (P)  storage  (Kornberg  1995;  Kornberg  et  al.  1999).  It  is  the  link  between 
polyP  and  phosphate  storage  that  has  driven  research  on  polyP  dynamics  as  a  function  of 
microbial  P  physiology. 

There  are  two  basic  processes  that  have  been  studied  involving  polyP  formation  in 
baeteria  in  response  to  P.  The  first,  termed  luxury  uptake,  is  the  storage  of  P  as  polyP 
when  P  is  consistently  in  excess  relative  to  other  macronutrients  like  nitrogen  (e.g., 
growth  is  not  P  eontrolled)  and  uptake  exeeeds  instantaneous  growth  demands.  Luxury 
uptake  has  been  extensively  studied  in  waste  water  treatment  applications  (Crocetti  et  al. 

2000) ,  and  in  marine  and  freshwater  phytoplankton  eulture  experiments  (e.g.,  Bertilsson 
et  al.  2003;  Diaz  et  al.  2008).  Luxury  uptake  could  thus  drive  the  aecumulation  of  polyP 
in  coastal  systems  or  areas  where  P  is  in  exeess  relative  to  nitrogen.  Seeond,  and 
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fundamentally  distinct  from  luxury  uptake,  is  the  overplus  response,  where  P  deplete  cells 
accumulate  polyP  in  response  to  short  term  pulses  in  P  supply  (Boiler  et  al.  1992; 
Jacobson  and  Halmann  1982).  The  overplus  response  is  not  well  studied  in  marine 
phytoplankton,  or  marine  systems  in  general,  but  it  has  been  hypothesized  that  overplus 
could  drive  the  cellular  accumulation  of  polyP  in  low  P  systems  where  phytoplankton  are 
P  deficient,  but  experience  short  term  variations  in  their  local  P  environment  (Karl  and 
Bjorkman  2001). 

Several  different  taxa  of  marine  phytoplankton  are  known  to  accumulate  polyP, 
including  Skeletonema,  Thalassiosim,  Synechocystis,  Nostoc,  Calothrix  (Diaz  et  al.  2008; 
Mateo  et  al.  2006)  and  Trichodesmium,  a  nitrogen-fixing  marine  cyanobacteria  (Romans 
et  al.  1994).  Putative  genes  involved  in  polyP  synthesis  and  degradation  also  appear  to 
be  present,  and  expressed,  in  a  number  of  both  eukaryotic  (Dyhrman  et  al.  2006b)  and 
prokaryotic  phytoplankton  (Gomez-Garcia  et  al.  2003;  Martiny  et  al.  2006;  Hewson  et  al. 
2009).  Taken  together,  these  observations  suggest  that  many  phytoplankton  produce  and 
degrade  polyP,  but  despite  its  potential  significance  to  P  physiology,  and  P  cycling,  there 
are  only  a  few  studies  that  have  assayed  polyP  accumulation  in  a  marine  setting  (Diaz  et 
al.  2008;  Romans  et  al.  1994).  Furthermore,  the  presence  of  phytoplankton  polyP  has 
never  been  examined  in  a  chronically  low  P  environment,  such  as  the  Sargasso  Sea. 

In  this  study,  Trichodesmium  near-surface  populations  from  the  Sargasso  Sea 
were  assayed  with  solid-state  ^^P  nuclear  magnetic  resonance  (NMR)  spectroscopy  to 
examine  cellular  allocation  of  P  as  polyP  in  an  oligotrophic  system.  Results  are 
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compared  to  culture  data,  and  co- varying  field  parameters,  to  examine  potential  drivers  of 
intracellular  polyP  dynamics. 

Methods 

Culture  Conditions 

Trichodesmium  cultures  were  grown  in  RMP  medium  with  Sargasso  seawater  as 
previously  described  (Webb  et  al.  2001).  Species  examined  include  Trichodesmium 
erythraeum  IMSIOI,  T.  tenue  Z-1,  T.  thiebautii  II-3,  and  T.  spiralis  KAT,  T.  spp.  H9-4 
obtained  from  the  culture  collection  of  Dr.  John  Waterbury  at  the  Woods  Hole 
Oceanographic  Institution.  All  cultures  were  unialgal  but  not  axenic.  Cultures  were 
grown  at  25°C  on  a  shaker  table  with  daily  cycles  that  consisted  of  1  hr  at  13.5  //E  m'  s'  , 
10  h  at  33.7  //E  m'^  s'\  1  h  at  13.5  //E  m'^  s'\  and  12  h  of  dark.  Chlorophyll  a 
fluorescence  was  monitored  using  an  aquaflor  fluorometer  (Turner  Designs).  T. 
erythraeum  IMSlOlwas  grown  with  (+P)  and  without  (-P)  16  pmol  E'^  phosphoric  acid 
added  to  the  media  and  samples  were  taken  for  particulate  carbon  (C),  nitrogen  (N)  and  P 

o  1 

as  well  as  for  solid-state  P  NMR  spectroscopy  (see  below)  throughout  the  experiments. 

Field  Samples 

Samples  were  collected  during  a  transect  of  the  Sargasso  Sea  in  May  2006  (Pig.  1). 
Trichodesmium  colonies  were  collected  from  the  near  surface  (roughly  within  the  top  20 
m)  using  a  hand  held  130  pm  net.  Colonies  were  not  sorted  by  morphology,  but  the  puff 
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morphology  was  typically  most  abundant.  Single  colonies  were  picked  into  0.2  //m- 
filtered  local  surface  seawater,  and  then  transferred  into  fresh  0.2  /rm-filtered  water  to 
reduce  contamination  of  closely  associated  organisms.  Colonies  were  then  separated  for 
further  analyses  as  indicated  below. 

NMR  spectroscopy 

Colonies  from  the  field  (approximately  100)  or  Trichodesmium  cultures  were  filtered 
onto  a  47mm  5  pm  polycarbonate  filter.  Filters  were  dried  at  65°C  prior  to  analysis, 
according  to  conditions  empirically  validated  in  other  work  (Benitez-Nelson  et  al.  2004). 

o  1 

Solid-state  P  NMR  spectra  were  recorded  on  a  Varian  Inova  500  spectrometer  operating 
at  202.489  MHz  using  a  Doty  Scientific  4  mm/XC  magic  angle  spinning  (MAS)  probe. 
Bloch  decays  of  50  ms  were  collected  with  a  200  ppm  window  after  30  degree  excitation 
pulses.  Two  pulse  phase  modulation  (Mcgeorge  et  al.  1999)  proton  dipolar  decoupling 
with  a  field  strength  of  45  kHz  was  applied  during  acquisition,  and  a  MAS  speed  of  10 
kHz  was  used.  Spectra  were  fit  with  5  Lorentzian  lines  to  determine  the  relative  ratio  of 
the  species  (polyP  =  -23  ppm).  The  percentage  of  total  P  as  polyP  was  determined  from 
the  signal  at  -23  ppm  (Fig.  2).  This  is  likely  to  be  a  minimum  measurement  given  the 
assumption  that  all  polyP  has  a  peak  signal  at  -23  ppm.  Results  are  reported  as  a  %polyP 
of  total  P  or  as  the  polyP:C.  The  later  value  was  determined  as  the  fraction  polyP  of  total 
P  multiplied  by  the  ratio  of  total  P:C  in  the  sample  to  account  for  changes  in  the  absolute 
amount  of  P  in  the  sample  due  to  variations  in  colony  size  or  P  quota.  The  %polyP 

o  1 

values  reported  have  an  error  of  +5%,  similar  to  solid-state  P  NMR  studies  in  other 
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systems  (Diaz  et  al.  2008). 


Alkaline  Phosphatase  Activity  (APA) 

For  Trichodesmium  APA  assays,  3-5  colonies  were  filtered  onto  a  5/^m  polycarbonate 
filter.  Samples  were  processed  as  described  elsewhere  (Dyhrman  and  Ruttenberg  2006). 
Briefly,  10  //mol  L'^  6,8-difluoro-4-methylumbelliferyl  phosphate  (di-MUF-P) 
(Invitrogen)  was  added  to  each  sample  in  a  petri  dish  with  artificial  seawater  containing 
no  P,  and  fluorescence  was  measured  on  a  Fluostar  Optima  plate  reader  (BMG  Labtech) 
every  5-20  minutes  for  5  time  points,  within  the  linear  range  of  the  assay.  Previous 
kinetics  experiments  found  the  10  //mol  L'^  substrate  concentration  to  be  saturating  (data 
not  shown).  Standard  fluorescence  curves  were  generated  for  each  assay  using  6,8- 
difluoro-7 -hydroxy -4-methylcoumarin  (di-MUF) . 

Chemical  analyses 

For  C,  N  and  total  P  measurements  20  Trichodesmium  colonies  from  the  field,  were 
collected  onto  pre-combusted  GF/F  filters  and  stored  dried  at  65°C  in  pre-combusted  foil. 
Each  filter  was  split  for  analysis  of  C,  N,  and  total  particulate  P.  C  and  N  concentrations 
were  measured  using  a  Perkin  Elmer  2400  CHN  Element  Autoanalyzer  with  no  acid 
fuming  (Karl  et  al.  1991).  Total  particulate  P  was  measured  as  described  in  Benitez- 
Nelson  et  al.  (2004).  Soluble  reactive  phosphorus  (SRP)  was  measured  according  to  the 
MAGIC-SRP  method  (Karl  and  Tien  1992),  with  a  detection  limit  of  0.5  nmol  E'\  SRP 
concentrations  measured  here  are  referred  to  as  the  dissolved  inorganic  phosphorus  (DIP) 
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concentration.  The  total  dissolved  P  (TDP)  was  processed  according  to  the  persulfate 
oxidation  method  (Valderrama  1981).  Dissolved  organic  P  (DOP)  concentrations  were 
calculated  as  the  difference  between  TDP  and  DIP.  With  the  methods  used  herein,  polyP 
would  be  detected  within  the  DOP  pool  (Monaghan  and  Ruttenberg  1999). 

Statistics 

Pair-wise  coefficients  of  correlation  and p  values  were  calculated  between  all  parameters; 
Trichodesmium  polyPrC,  Trichodesmium  APA,  N:P  ratio,  C:P  ratio,  surface  temperature, 
mixed  layer  depth,  latitude,  DOP  and  DIP  concentration  at  20  m  using  IMP  software 
(SAS,  Cary  NC).  The  mixed  layer  depth  was  defined  as  where  the  temperature  changed 
by  0.5°C  relative  to  the  surface.  A  two-tailed  t-test  was  used  to  compare  %polyP  among 
growth  phases  of  cultured  Trichodesmium  and  between  -i-P  and  -P  culture  treatments. 

Genome  annotation 

The  putative  ppX  and  ppK  genes  were  identified  using  the  Integrated  Microbial  Genomes 
(IMG)  portal  at  the  Joint  Genome  Institute-  US  Department  of  Energy  website 
(http://img.jgi.doe.gov/vI .  I/main.cgi)  and  the  National  Center  for  Biotechnology 
Information  Genbank  database  (http://www.ncbi.nlm.nih.gov/),  based  on  homology. 

DNA  Amplification  and  Sequencing 

Cultures  of  Trichodesmium  were  collected  by  filtration  onto  a  5  pm  polycarbonate  filter 
and  stored  at  -20°C  until  extraction.  DNA  was  extracted  using  the  Instagene  Matrix  (Bio- 
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Rad)  according  to  the  manufacturer’s  instructions  or  as  described  in  Ehrenreich  et  al. 
(2005).  PCR  primers  were  designed  to  either  amplify  fragments  of  the  putative  ppX  (F  5’ 
GGAATGTCCGAAAAGCGAGC  3’  R  5’  GCCCAAAAAGCAACCCCGTTC  3’)  and 
ppK  (F  5’  CGCTTTATCAAACTGATTCGTCGG  3’ 

R  5’  CGCAAACAACAAATACCACGGAC  3’)  genes.  Each  PCR  reaction  consisted  of 
a  5  min  denaturation  step  at  95°C,  followed  by  35  cycles  of  1  min  at  95°C,  1  min  at 
60.7°C,  1  min  at  72°C  and  a  final  extension  of  10  min  at  72°C,  in  a  Bio-Rad  iCycler. 
PCR  amplification  was  done  using  0.5  //F  iTaq  DNA  polymerase  (Bio-Rad,  Hercules, 
CA),  2.5  /tF  of  DNA  template,  2.5/tF  of  2  mM  dNTP,  2.5  /tF  of  5x  iTaq  buffer,  50  pmol 
of  each  primer,  and  sterile  water  to  a  final  volume  of  25  /tF.  PCR  products  were  gel 
extracted  using  the  QIAquick  gel  extraction  kit  (Qiagen),  and  direct  sequenced  at  the 
MWG  Biotech  facility  or  at  the  University  of  Maine,  according  to  the  facility’s  protocols. 
Sequences  were  edited  using  Sequencher  (Gene  Codes  Inc.)  and  verified  manually. 
Alignments  were  done  in  Mac  Vector.  The  sequences  have  been  deposited  in  Genbank 
(Accession#  GU299287-GU299292). 

Results 

Genes 

T.  erythraeum  has  gene  homologs  for  a  ppK,  and  for  ppX.  The  genes  are  not  contiguous, 
and  do  not  appear  to  be  downstream  of  a  Pho  box  P  regulatory  sequence  (Su  et  al.  2007), 
although  previous  work  has  highlighted  that  there  may  be  heterogeneity  in  the 
Trichodesmium  Pho  box  sequence  (Orchard  et  al.  2009).  All  species  of  Trichodesmium 
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tested  {T.  erythraeum,  T.  spiralis,  T.  tenue,  and  T.  thiebautii)  had  ppK  and  ppX.  For  these 
species  over  the  sequenced  156  base  pair  fragment,  ppK  was  88-93%  identical  at  the 
nucleotide  level,  and  the  sequenced  223  base  pair  fragment  of  ppX  was  93-100%  identical 
at  the  nucleotide  level. 

Culture  experiments 

Multiple  Trichodesmium  species  had  detectable  polyP.  The  maximum  %polyP  that  was 
measured  in  P  replete  cultures  was  12%  for  T.  tenue,  14%  for  T.  spp.  H9-4,  19%  for  T. 
thiebauttii  and  22%  for  T.  erythraeum.  The  %polyP  in  P  replete  T.  erythraeum  IMSIOI 
cultures  was  significantly  lower  in  log  phase,  than  in  stationary  phase  (Table  1).  The 
%polyP  in  log  phase  P  replete  cultures  (cultures  with  a  C:P  ratio  less  than  200)  was 
significantly  lower  than  log  phase  P  deplete  cultures  (grown  with  no  added  P  and  with  a 
C:P  ratio  greater  than  200)  (Table  1). 

Field  data 

The  surface  temperature  increased  with  decreasing  latitude  from  north  to  south  increasing 
from  24.6°C  to  27.7°C  (Table  2).  The  mixed  layer  depth  generally  shoaled  in  the  more 
southern  latitudes  decreasing  from  a  maximum  of  32  m  to  a  minimum  of  13  m  (Table  2). 
DIP  and  DOP  concentrations  at  20  m  were  relatively  constant  along  the  cruise  transect, 
ranging  from  0.5  nmol  L'*-  3.0  nmol  L'^  and  22.9  nmol  L'^-  46.8  nmol  L'^  respectively, 
and  with  no  clear  trend  from  north  to  south  (Table  2).  The  average  Trichodesmium  C:P 
ratio  and  N:P  ratio  were  above  the  Redfield  ratio  of  106:1  and  16:1  respectively  and 
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Trichodesmium  APA  ranged  from  0.23  -  0.97  nmol  P  hr'^  colony'\  None  of  these 
potential  metrics  of  P  physiology  showed  a  consistent  trend  on  the  north  to  south  transect 
(Table  2),  or  a  significant  correlation  to  other  environmental  parameters.  However,  these 
metrics  of  P  physiology  were  consistently  in  the  range  of  what  would  be  considered 
evidence  of  P  depletion  from  other  studies  (Krauk  et  al.  2006;  White  et  al.  2006^?). 

In  field-collected  Trichodesmium  colonies,  %polyP  ranged  from  8  -  25%  of  the 

T  1 

total  P  as  detected  with  solid-state  P  NMR  spectroscopy  (Table  2).  The  polyP 
normalized  between  stations  as  the  polyP:C  ratio,  was  significantly  correlated  with 
temperature,  latitude,  and  mixed  layer  depth  (Table  2).  There  was  no  significant 
correlation  between  Trichodesmium  polyP:C  and  DOP,  DIP  concentration, 
Trichodesmium  N:P  ratio,  or  Trichodesmium  APA.  There  was  also  no  significant 
correlation  between  APA  and  DIP  concentration,  C:P  ratio  or  N:P  ratio. 

Discussion 

Storage  of  cellular  P  as  polyP,  is  considered  an  important  aspect  of  microbial  P 
physiology,  and  could  influence  P  cycling  in  aquatic  systems.  However,  there  are  very 
few  studies  that  have  examined  polyP  in  marine  phytoplankton,  and  none  have  examined 
polyP  accumulation  in  low  P  oligotrophic  systems.  Others  have  hypothesized  that 
phytoplankton  polyP  could  accumulate  in  phytoplankton  from  oligotrophic  low  P  regions 
as  part  of  an  overplus  type  response  (Karl  and  Bjorkman  2001).  To  address  this 
hypothesis  we  examined  polyP  in  Trichodesmium  cultures  and  in  field  populations  from 
the  Sargasso  Sea. 
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PolyP  dynamics  in  culture 

All  of  the  Trichodesmium  species  tested  in  culture  have  the  genetic  machinery  to 

o  1 

produce  and  degrade  polyP.  This  is  consistent  with  solid-state  P  NMR  data  that 
detected  polyP  in  P  replete  cultures  of  Trichodesmium,  with  maximal  polyP  percentages 
that  were  similar  among  species.  These  data  suggest  that  the  polyP  dynamics  observed  in 
T.  erythraeum  IMS  101  are  likely  representative  of  the  other  species,  although  this 

o  1 

warrants  further  study.  The  solid-state  P  NMR  method  is  non-destructive  and  can 
provide  data  on  P  allocation  patterns  that  is  complimentary  to  other  approaches  that  have 
detected  polyP  accumulated  in  granules  in  Trichodesmium  populations  from  the 
Caribbean  Sea  using  electron  microscopy  (Romans  et  al.  1994). 

Under  P  replete  conditions  T.  erythraeum  IMS  101  had  the  highest  %polyP  of  total 
cellular  P  during  stationary  phase.  In  fact,  polyP  was  often  below  the  detection  limit  in  P 
replete  log  phase  cultures.  This  is  consistent  with  the  hypothesis  that  polyP  is  important 
for  survival  in  stationary  phase  or  that  polyP  is  accumulating  in  response  to  general  stress 
(Rao  and  Kornberg  1996;  Kornberg  and  Ault-Riche  1999).  Luxury  uptake  may  also  be 
occurring  in  these  cultures,  but  luxury  P  is  either  not  being  stored  as  polyP  or  is  below  the 
detection  limit  for  this  assay.  These  results  contrast  with  those  from  P  deplete  cultures, 
where  log  phase  P  deplete  samples  had  significantly  higher  percentages  of  polyP  relative 
to  log  phase  P  replete  samples.  This  could  be  the  result  of  an  overplus-like  response. 
Because  DIP  was  not  re-fed  to  P  deplete  cells  this  is  not  entirely  consistent  with 
traditional  definitions  of  the  overplus  response.  However,  upon  becoming  P  deplete  T. 
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erythraeum  up-regulates  a  suite  of  genes  encoding  enzymes  to  access  P  from  DOP 
(Dyhrman  et  al.  2006a;  Orchard  et  al.  2009).  These  cultures  were  grown  on  Sargasso  Sea 
water  containing  DOP,  and  cells  would  have  been  able  to  access  a  new  source  of  P  upon 
DIP  depletion.  This  DOP  access  and  utilization  might  have  resulted  in  an  overplus  like 
accumulation  of  polyP.  This  could  be  analogous  to  what  could  happen  in  P  deplete 
oligotrophic  systems  where  cells  experience  variation  in  the  P  supply  while  switching 
between  growth  on  low  concentration  DIP,  to  growth  on  comparatively  higher 
concentrations  of  DOP. 

PolyP  dynamics  in  the  Sargasso  Sea 

PolyP  accumulation  in  phytoplankton  has  previously  been  observed  in  the  coastal 
zone  and  eutrophic  lakes,  and  in  these  regions  polyP  accumulation  was  generally 
attributed  to  luxury  uptake  of  P  (Diaz  et  al.  2008;  Schelske  and  Sicko-Goad  1990).  In  the 
oligotrophic  Sargasso  Sea  polyP  was  between  8-25%  of  the  total  cellular  P  in 
Trichodesmium  despite  the  low  DIP  concentrations  in  this  area  (<5  nmol  L'*).  These  data 
are  striking  in  that  they  demonstrate  that  even  in  low  P  regions  some  phytoplankton 
sequester  a  significant  proportion  of  their  P  as  polyP.  These  field  data  are  consistent  with 
the  culture  results  that  highlight  the  allocation  of  P  to  polyP  can  be  elevated  in  P  deplete 
T richodesmium . 

Unlike  the  other  marine  systems  that  have  been  examined  (e.g.  Carribean; 
Romans  e  al.  1994),  the  presence  of  high  percentages  of  polyP  in  Trichodesmium  from 
the  Sargasso  Sea  is  unlikely  to  be  due  to  luxury  uptake.  Evidence  from  culture  studies 
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suggest  Trichodesmium  may  be  capable  of  luxury  uptake  (White  et  al.  2006^?),  but 
elevated  N:P,  C:P,  and  APA  (metrics  of  P  deficiency)  in  this  study  and  others  from  the 
region  suggest  that  Trichodesmium  populations  are  P  deficient  (Dyhrman  et  al.  2002; 
Mulholland  et  al.  2002;  Sohm  et  al.  2008).  Luxury  uptake  is  probably  not  a  driving 
factor  for  polyP  accumulation  in  either  Trichodesmium  or  other  P  deficient  phytoplankton 
from  low  P  oligotrophic  regimes. 

Trichodesmium  is  unique  relative  to  some  other  phytoplankton  in  that  some 
species  can  regulate  buoyancy.  It  has  been  hypothesized  that  this  organism  may  migrate 
to  the  phosphocline,  store  P  as  polyP,  and  return  to  the  surface  in  a  process  termed  P- 
mining  (Villareal  and  Carpenter  2003).  Although  this  could  result  in  elevated  polyP,  the 
phosphocline  ranged  from  160-300  m  on  this  transect  which  is  deeper  than 
Trichodesmium  is  likely  to  be  able  to  migrate  (White  et  al.  2006a),  and  deeper  than 
Trichodesmium  was  observed  in  this  region  (Davis  and  McGillicuddy  2006;  Orcutt  et  al. 
2001).  Additionally,  a  study  from  the  Caribbean  found  that  polyP  granules  in 
Trichodesmium  were  more  abundant  in  negatively  buoyant  colonies  than  in  positively 
buoyant  colonies;  a  trend  opposite  of  that  expected  if  P-mining  were  occurring  (Romans 
et  al.  1994).  These  factors  suggest  that  P-mining  is  not  driving  the  high  percentages  of 
polyP  observed  in  Trichodesmium  from  the  Sargasso  Sea. 

In  the  field,  Trichodesmium  may  increase  allocation  of  P  to  polyP  due  to  an 
overplus-like  response  upon  up-regulating  genes  involved  in  DOP  acquisition,  or 
experiencing  a  change  in  P  supply  in  an  already  P  deficient  population.  There  is  the 
additional  possibility  that  in  the  environment  Trichodesmium  colonies  are  encountering 
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micro-scale  variability  in  P  availability,  and  thus  constantly  experiencing  fluxes  in  P 
supply.  These  micro-scale  fluctuations  in  P  at  the  single  cell  or  colony  level  are  not 
captured  by  traditional  bulk  chemical  sampling.  The  importance  of  understanding  the 
environment  of  the  microorganism  at  the  nanometer  or  micrometer  scale  is  increasingly 
appreciated  (Azam  and  Malfatti  2007),  and  is  likely  to  be  crucial  to  understanding  the 
nutritional  physiology  of  phytoplankton  in  situ. 

Comparisons  between  PolyP  and  other  environmental  factors 

To  examine  how  physical,  chemical  and  physiological  parameters  influence  polyP 
formation  along  the  transect,  Trichodesmium  polyP:C  was  calculated  to  account  for 
changes  in  the  absolute  quantity  of  polyP  and  total  P.  PolyP:C  was  not  correlated  to 
Trichodesmium  physiology  as  indicated  by  C:P,  N:P  or  APA  activity,  although,  with  the 
possible  exception  of  Sta.  3,  all  samples  were  clearly  P  deplete  based  on  these  indicators. 
No  gradient  in  P  physiology  was  observed  across  stations  that  had  such  low  DIP 
concentrations.  PolyP:C  was  also  not  correlated  to  DOP  or  DIP  concentration,  and 
therefore  polyP:C  does  not  seem  to  be  dependant  on  the  standing  stock  concentration  of 
P.  However,  the  DIP  and  DOP  concentrations  were  very  low  across  the  transect,  and  it 
may  be  that  comparisons  with  other  systems  with  higher  DIP  concentrations  are 
necessary  to  resolve  differences  in  P  allocation  as  a  function  of  P  standing  stock.  APA  is 
known  to  be  correlated  with  large  changes  in  DIP  concentration  in  previous  studies 
(Mulholland  et  al.  2002;  Dyhrman  and  Ruttenberg  2006).  However,  in  this  study  it  was 
not,  further  underscoring  the  fact  that  there  was  not  a  gradient  in  P  concentrations  and 
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Trichodesmium  physiology  because  of  the  consistently  low  DIP  concentrations.  The  lack 
of  relationship  is  consistent  with  the  hypothesis  that  micro-scale  flux  in  P,  or  shifts 
between  DIP  and  DOP  metabolism,  in  a  P  deplete  environment,  is  driving  polyP 
accumulation,  not  the  overall  standing  stock. 

There  was  a  correlation  between  polyP:C  and  latitude  and  other  physical 
parameters  that  co-vary  with  latitude  such  as  temperature  and  mixed  layer  depth,  where 
polyP:C  was  the  highest  at  the  northern  end  of  the  transect  with  cooler  temperatures,  and 
a  deeper  mixed  layer.  Work  with  other  microbes  found  that  polyP  concentrations 
increased  with  heat  stress,  or  with  general  stress  (Kornberg  et  al.  1999;  Seufferheld  et  al. 
2008).  In  this  study,  polyP  accumulated  with  decreasing  temperature,  which  is  the 
opposite  of  the  previously  observed  pattern,  but  may  indicate  increased  stress  in  the  more 
northern  stations,  given  that  Trichodesmium" s  growth  is  optimum  at  24-30°C  (Breitbarth 
et  al.  2006).  Therefore,  the  northern  most  stations  were  at  the  edge  of  the  optimal  growth 
temperatures  for  Trichodesmium.  Clearly  many  different  factors  can  influence  polyP  and 
it  would  be  useful  to  examine  polyP  in  Trichodesmium  cultures  and  field  populations 
under  a  greater  suite  of  conditions  to  determine  the  consistency  of  this  trend  in  the  field 
data. 

Impact  of  polyP  to  the  biogeochemistry  of  the  Sargasso  Sea 

It  remains  to  be  determined  if  elevated  P  allocation  to  polyP  is  a  trait  common  to 
phytoplankton  in  regions  like  the  Sargasso  Sea,  or  unique  to  Trichodesmium,  yet  this  can 
have  profound  ramifications  for  the  P  cycle  of  oligotrophic  systems.  It  is  widely 
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accepted  that  microbial  phosphoesters  (Ammerman  et  al.  2003;  Benitez-Nelson  2000; 
Benitez-Nelson  and  Buesseler  1999),  and  potentially  microbial  phosphonates  (Dyhrman 
et  al.  2009)  contribute  to  upper  water  column  DOP  and  microbial  polyP  may  also  be 
present  and  dissolved  in  the  upper  water  column.  While  polyP  is  an  inorganic 
compound,  known  to  be  bioavailable  to  support  production  (Moore  et  al.  2005),  it  is  often 
measured  as  DOP,  and  could  represent  a  previously  ignored  P  source  for  production  in 
the  upper  water  column.  Diaz  et  al  (2008)  observed  a  similar  %polyP  in  phytoplankton 
and  in  DOP  in  a  coastal  region  (7%  and  11%  respectively).  Assuming  Trichodesmium 
polyP  accumulation  is  representative  of  Sargasso  Sea  phytoplankton  in  general,  then  we 
can  extrapolate  from  Trichodesmium  to  the  water  column.  On  average  in  the  Sargasso 
Sea  there  was  15%  polyP  in  Trichodesmium.  If  the  average  of  49  nmol  L'^  DOP  for  this 
transect  is  15%  polyP,  then  the  standing  stock  of  polyP  in  the  upper  water  column  is  on 
average  5  nmol  L'\  This  is  higher  than  the  average  DIP  concentration  of  1.2  nmol  L"', 
and  would  more  than  double  the  bioavailable  inorganic  P  supply.  Further  study  is 
warranted  to  assay  for  the  presence  of  polyP  in  phytoplankton  and  the  water  column  of 
this  region,  given  the  possible  importance  of  polyP  to  the  P  cycle  in  this  and  potentially 
other  oligotrophic  systems. 

Conclusions 

In  this  study,  we  found  that  Trichodesmium  forms  polyP  in  a  low  P  environment 
(the  Sargasso  Sea),  and  that  polyP  can  be  a  major  component  of  Trichodesmium  P  pool. 
While  it  is  unclear  what  exact  factors  induce  the  formation  of  polyP  in  this  system,  the 
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data  are  consistent  with  a  possible  overplus-type  response.  It  may  be  that  Trichodesmium 
is  taking  advantage  of  fluxes  in  P  supply,  and  that  polyP  is  a  physiological  fingerprint  of 
this  variability.  Additionally,  if  polyP  formation  is  a  common  trait  in  phytoplankton, 
these  data  suggest  that  polyP  might  be  an  important  additional  bio-available  component 
of  DOP.  Taken  together,  this  study  highlights  the  potential  importance  of  this  compound, 
and  suggests  the  importance  of  further  studies  into  the  dynamics  of  polyP  formation  in 
oligotrophic  systems. 
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Table  1.  The  %polyP  detected  in  T.  erythraeum  IMS  101  cultures  grown  on  phosphorus 


(P)  replete  or  P  deplete  media  in  different  growth  phases. 


Log  phase 

P  replete^ 

(n=8) 

Stationary  phase 

P  replete 
(n=4) 

Log  phase 

P  deplete 
(n=3) 

mean 

0.2’’ 

16 

16 

standard  deviation 

0.5 

9 

13 

p  value‘s 

0.0003 

0.006 

‘T'he  average  C:P  ratio  for  the  P  replete  cultures  was  130+35,  and  was  486+201  for  the  P 
deplete  samples. 

'’Each  %polyP  value  has  a  +5%  error. 

‘^All  samples  were  compared  to  the  log  phase  P  replete  condition  using  a  two-tailed  t  test, 
bold  font  denotes  significance. 
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Each  %polyP  value  has  a  +5%  error. 

2 

Pairwise  correlations  are  compared  to  polyP:C  at  each  station,  bold  font  denotes  significance 
nd  is  not  determined. 
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Figure  1.  Station  locations  for  a  May  2006  transect  from  north  to  south  in  the  Sargasso 
Sea. 
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Figure  2.  Representative  solid-state  NMR  spectrum  from  Station  2  during  the  May 
2006  cruise  transect  in  the  Sargasso  Sea,  the  polyphosphate  peak  is  noted  at  -23  ppm. 
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Summary 

The  marine  diazotroph  Trichodesmium  is  a  major 
contributor  to  primary  production  and  nitrogen  fixa¬ 
tion  in  the  tropicai  and  subtropicai  oceans.  These 
regions  are  often  characterized  by  low  phosphorus 
(P)  concentrations,  and  P  starvation  of  Trichodes¬ 
mium  couid  limit  growth,  and  potentiaiiy  constrain 
nitrogen  fixation.  To  better  understand  how  this 
genus  responds  to  P  starvation  we  examined  four 
genes  invoived  in  P  acquisition:  two  copies  of  a  high- 
affinity  phosphate  binding  protein  (pstS  and  sphX) 
and  two  putative  aikaiine  phosphatases  {phoA  and 
phoX).  Sequence  analysis  of  these  genes  among  cui- 
tured  species  of  Trichodesmium  (T.  tenue,  T.  eryth- 
raeum,  T.  thiebautii  and  T.  spiraiis)  showed  that  they 
ali  are  present  and  conserved  within  the  genus.  In 
T.  erythraeum  IMS101,  the  expression  of  sphX,  phoA 
and  phoX  were  sensitive  to  P  suppiy  whereas  pstS 
was  not.  The  induction  of  aikaiine  phosphatase  activ¬ 
ity  corresponded  with  phoA  and  phoX  expression, 
but  enzyme  activity  persisted  after  the  expression  of 
these  genes  returned  to  basai  leveis.  Additionaliy, 
/7/yH(nitrogenase  reductase;  involved  in  nitrogen  fixa¬ 
tion)  expression  was  downreguiated  under  P  starva¬ 
tion  conditions.  These  data  highlight  molecular  level 
responses  to  low  P  and  lay  a  foundation  for  better 
understanding  the  dynamics  of  Trichodesmium  P 
physiology  in  low-P  environments. 


Introduction 

Photosynthetic  phytoplankton  plays  a  key  role  in  the  bio¬ 
geochemistry  of  the  ocean  by  fixing  CO2  and  influencing 
the  flux  of  carbon  to  the  deep  sea.  The  export  of  this 
carbon  impacts  the  atmospheric  CO2  concentration,  and 
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therefore  climate,  over  geological  time  scales  (Longhurst, 
1991).  Photosynthetic  diazotrophs,  such  as  the  cyano¬ 
bacterium  Trichodesmium,  are  of  particular  importance  to 
the  biogeochemistry  of  the  ocean  because  they  fix  both 
carbon  and  nitrogen  (N).  Thus,  they  contribute  to  carbon 
cycling  directly,  as  primary  producers,  and  indirectly  by 
fixing  atmospheric  N,  which  can  then  promote  the  growth 
of  other  phytoplankton  species  (Capone,  2001 ;  Coles 
etal.,  2004).  For  example,  the  input  of  ‘new  N’  fixed  by 
Trichodesmium  has  been  estimated  to  be  as  much  as 
1.6  X  10^^  mol  N  into  the  tropical  North  Atlantic  (Capone 
etal.,  2005). 

The  factors  that  constrain  Trichodesmium  growth  and 
N2  fixation  are  intensely  studied  due  to  the  aforemen¬ 
tioned  importance  of  this  genus  to  global  biogeochemistry. 
Although  iron  is  often  thought  to  limit  N2  fixation  by 
Trichodesmium  (Berman-Frank  2001;  Webb  etal.,  2001; 
Kustka  etal.,  2002;  Mills  etal.,  2004),  recent  studies  in 
the  Atlantic  suggest  that  phosphorus  (P)  may  also  be 
limiting  (Sahudo-Wilhelmy  etal.,  2001;  Dyhrman  etal., 
2002;  Mulholland  etal.,  2002;  Webb  etal.,  2007).  Further 
supporting  these  observations,  dissolved  inorganic  phos¬ 
phate  (DIP)  concentrations  in  the  Atlantic  are  significantly 
lower  than  those  found  in  a  number  of  other  oligotrophic 
regimes,  such  as  the  Pacific  (Wu  etai,  2000;  Cavendar- 
Bares  etal.,  2001). 

The  major  forms  of  P  in  the  marine  environment  are  DIP 
and  dissolved  organic  phosphorus  (DOP).  DIP  is  consid¬ 
ered  the  preferable  P  source  for  bacteria  and  phytoplank¬ 
ton  because  it  is  directly  available  for  growth  (Vershinina 
and  Znamenskaia,  2002;  Dyhrman  et  a!.,  2007),  but  many 
microbes  are  also  able  to  acquire  P  from  DOP  and  this 
ability  has  been  observed  in  a  number  of  marine  cyano¬ 
bacteria  and  other  phytoplankton  (Dyhrman  etal.,  2002; 
Dyhrman  and  Palenik,  2003;  Moore  etal.,  2005).  In  the 
upper  water  column  of  oligotrophic  regimes,  DOP  often 
comprises  a  significant  portion  of  the  total  P  (Orrett  and 
Karl,  1987;  Karl,  1999;  Wu  etal.,  2000),  and  a  recent 
study  by  Mather  and  colleagues  (2008)  suggests  that 
DOP  could  fuel  up  to  30%  of  the  primary  productivity  in  the 
North  Atlantic  (Mather  etal.,  2008).  Within  the  high- 
molecular-weight  DOP  pool  that  is  amenable  to  charac¬ 
terization  (approximately  20-40%),  roughly  25%  has  a 
phosphonate  bond  and  roughly  75%  has  an  ester  bond 
(Clark  etal.,  1998).  Therefore,  P  esters  represent  a  larger 
potential  P  source  than  phosphonates.  It  has  been  shown 
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that  Trichodesmium  is  capable  of  growing  on  model  phos- 
phomonoesters  (Stihl  etai,  2001),  and  gene  expression 
analysis  has  shown  that  phosphonate  compounds  are 
also  likely  available  to  this  genus  (Dyhrman  etal.,  2006). 
While  the  genes  for  phosphonate  acquisition  have  been 
examined,  the  genes  encoding  the  alkaline  phosphatase 
enzyme  responsible  for  the  hydrolysis  of  the  larger,  P 
ester  pool  have  not  been  studied  in  Trichodesmium. 

In  order  to  survive  in  low-P  environments,  microbes 
typically  induce  a  set  of  genes  involved  in  high-affinity 
phosphate  uptake  and  DOP  hydrolysis.  P  uptake  in  many 
microbes  occurs  via  the  P-specific  transport  (Pst)  system 
during  P-limiting  conditions,  or  via  the  P-inorganic  trans¬ 
port  (Pit)  system  under  P-replete  conditions  (Vershinina 
and  Znamenskaia,  2002).  DOP  hydrolysis  can  be 
mediated  by  many  different  enzymes,  including  alkaline 
phosphatases,  encoded  by  the  phoX  and  phoA  genes 
(Roy  etal.,  1982;  Moore  etal.,  2005;  Wu  etal.,  2007). 
In  Escherichia  coli  and  Synechocystis  PCC6803  these 
genes  are  under  the  control  of  a  histidine-kinase  tran¬ 
scriptional  regulator  (Vershinina  and  Znamenskaia,  2002; 
Su  etal.,  2003).  This  response  regulator  binds  to  a 
sequence  referred  to  as  the  Pho  box  in  the  promoter 
region  of  each  target  gene  or  gene  cluster.  Searching  for 
Pho  boxes  has  proved  to  be  a  useful  approach  for  iden¬ 
tifying  genes  that  may  be  involved  in  P  acquisition  (Su 
et  al.,  2007),  although  the  regulation  and  function  of  these 
genes  must  be  empirically  validated. 

Trichodesmium  spp.  are  dominant  members  of  the  phy¬ 
toplankton  community  in  tropical  and  subtropical  oceans 
(Capone  et  al.,  1 997),  with  six  described  species;  T.  eryth- 
raeum,  T.  tenue,  T.  thiebautii,  T.  hildebrandtii,  T.  contor- 
tum  and  T.  spiralis,  formerly  known  as  Katagnymene 
spiralis  (Janson  etal.,  1995;  Lundgren  etal.,  2005). 
These  species  differ  both  in  morphology  and  distribution, 
with  T.  thiebautii  and  T.  erythraeum  being  the  species 
most  frequently  associated  with  blooms  in  the  open  ocean 
(Orcutt  etal.,  2002).  In  this  study,  we  examine  T.  eryth¬ 
raeum,  T.  tenue,  T.  thiebautii  and  T.  spiralis  for  the  pres¬ 
ence  of  two  genes  involved  in  DIP  uptake  (the  putative  P 
binding  proteins:  pstS  and  sphX),  and  two  genes  involved 
in  DOP  hydrolysis  (the  putative  alkaline  phosphatases: 
phoA  and  phoX).  Using  T.  erythraeum  as  a  model,  we 
also  examine  how  these  genes,  and  one  gene  involved  in 
N2  fixation  (Nitrogenase  reductase:  nifH),  are  regulated  in 
response  to  changes  in  P  supply. 


Results 

Genomio  observations 

The  genome  of  T.  erythraeum  IMS101  contains  homo- 
logues  for  genes  involved  in  the  acquisition  of  DIP,  and  the 
hydrolysis  of  DOP  (Table  1).  Homologues  to  each  of  the 


genes  found  in  the  Pst  cluster  (PstCAB)  and  two  homo¬ 
logues  of  the  high-affinity  phosphate  binding  protein,  pstS 
(YP_723097)  and  sphX  (YP_723095)  are  present  in  the 
T.  erythraeum  genome  (Table  1).  These  two  genes  are 
separated  by  a  small  gap,  and  the  sphX  gene  is  preceded 
by  a  putative  Pho  box  sequence,  whereas  the  pstS  gene 
is  not  (Su  etal.,  2007)  (Fig.  1).  There  are  no  clear  homo¬ 
logues  to  any  genes  in  the  Pit  system. 

There  are  three  putative  alkaline  phosphatases  in 
the  T.  erythraeum  genome.  One  putative  alkaline  phos¬ 
phatase,  phoA  (YP_723031),  is  related  to  an  atypical 
alkaline  phosphatase  in  Synechococcus  (Ray  etal., 
1991),  to  which  it  has  a  52%  identity  over  approximately 
50%  of  the  translated  gene.  The  other  two  genes,  phoX 
(YP_723360)  and  phoX2  (YP_723924),  which  encode 
putative  alkaline  phosphatases,  share  the  same  con¬ 
served  protein  domain  (COG  3911)  as  the  characterized 
PhoX  proteins  of  Pseudomonas,  Pasterella  and  Vibrio 
(Pho''°  in  Vibrio)  (Roy  etal.,  1982;  Monds  etal.,  2006;  Wu 
etal.,  2007).  Both  phoX  and  phoX2  have  a  39%  identity 
over  99%  of  the  translated  gene  to  Pho'^'^  in  Vibrio.  A  Pho 
box  has  been  identified  in  front  of  the  putative  phoX  gene 
(Su  etal.,  2007),  but  not  before  the  putative  phoA  or 
phoX2  genes  (Table  1 ).  The  phoA  gene  does  have  a  very 
similar  repeat  to  that  of  a  Pho  box  in  its  promoter  region, 
but  lacks  the  three  base  linker  sequence  common  to  the 
Pho  boxes  identified  in  Su  and  colleagues  (2007) 
(Table  1).  Further  analyses  focused  on  phoA  and  phoX. 

Sequence  analysis 

The  genes  pstS,  sphX,  phoX  and  phoA  are  present  and 
highly  conserved  in  all  Trichodesmium  species  examined. 
With  the  exception  of  pstS,  all  genes  were  96-100% 
identical  at  the  nucleotide  level  over  the  gene  fragments 
sequenced  (Table  2).  For  the  pstS  gene  there  appear  to 
be  two  groups:  T.  tenue  and  T.  erythraeum,  which  are 
99%  identical,  and  T.  spiralis  and  T.  thiebautii,  which  are 
98%  identical,  but  the  two  pstS  groups  are  only  87-88% 
identical  to  each  other  (Table  2). 

Culture  experiments 

The  expression  of  genes  involved  in  P  uptake  and  DOP 
hydrolysis  were  examined  in  two  ways:  (i)  by  transferring 
exponentially  growing  cells  to  replete  (-rP)  and  no  P  (-P) 
media,  in  the  transfer  experiment  or  (ii)  by  starving  cells 
for  P  and  then  adding  back  P,  in  the  re-feed  experiment. 
Each  of  these  experiments  was  replicated  with  similar 
results.  Here  we  present  the  data  from  one  of  each  of 
these  experiments. 

In  the  transfer  experiment,  alkaline  phosphatase  activity 
(APA)  increased  in  the  -P  culture  to  a  final  rate  of 
16.5  nmol  P  h“^  pg  Chi  a~\  and  was  statistically  greater 
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P-regulated  genes  were  selected  based  on  available  regulation  information,  and  relevance  to  the  Trichodesmium  genes  examined  in  this  study.  Several  genes  are  listed  that  are  part  of  the  same 
gene  cluster,  and  therefore  are  only  preceded  by  one  Pho  box  for  both  genes.  In  those  instances  the  Pho  box  sequence  was  only  listed  for  the  first  gene  in  the  gene  cluster. 

NE;  not  examined. 
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Fig.1.  Pst  topology  of  Trichodesmium  erythraeum  IMS101  and  Synechocystis  PCC6803  with  the  percent  identity  reported  for  amino  acid 
alignments  over  80%  of  the  translated  gene.  The  pstS  and  sphX  genes  encode  high-affinity  binding  proteins,  pstA,  pstB  and  pstC  are 
components  of  the  phosphate  transporter.  A  vertical  arrow  next  to  the  gene  indicates  the  location  of  a  Pho  box  sequence  as  identified  in  Su 
and  colleagues  (2007).  Genes  shaded  grey  have  been  shown  to  be  upregulated  by  P  starvation  in  this  study  or  in  Suzuki  and  colleagues 
(2004). 


than  the  -hP  cultures  at  168  h  (P<  0.05)  (Fig.  2).  In  the  -hP 
cultures,  APA  reached  a  maximum  level  of  4.0  nmol 
P  h"^  |4g  Chi  a~\  The  expression  of  phoA  and  phoX 
tracked  with  the  measured  APA.  The  expression  of  phoA 
in  the  -P  culture  was  statistically  different  from  the  -hP 
culture  by  72  h  (P=  0.031),  as  was  the  expression  of 
phoX  (P=  0.028)  (Fig.  2).  In  the  re-feed  experiment,  APA 
started  at  24.9  nmol  P  h“^  pg  Chi  a~',  this  rate  declined 
upon  P  addition,  and  was  statistically  different  from  the  -P 
culture  by  48  h  (P<  0.05)  (Fig.  2).  The  phoA  and  phoX 
genes  were  downregulated  compared  with  the  -P  culture 
at  each  time  point  (Fig.  2).  This  difference  was  statistically 
significant  by  4  h  for  phoX  (P<  0.001)  and  phoA 
(P<  0.001)  (Fig.  2). 

The  expression  of  sphX  was  upregulated  in  the  -P 
cultures  in  the  transfer  experiment,  and  was  statistically 
different  from  the  -hP  culture  at  the  216  h  time  point 
(P=  0.021).  While  pstS  expression  was  always  detect¬ 
able,  there  was  no  significant  change  in  expression  during 
the  time-course  of  this  experiment  (Fig.  3).  In  the  re-feed 
experiment,  sphX  was  downregulated  in  response  to  the 
P  addition  (Fig.  3).  The  difference  in  expression  between 
the  -P  and  re-fed  cultures  was  statistically  significant  by 
24  h  (P=  0.036).  There  was  no  clear  pattern  in  expression 


for  pstS,  but  it  was  significantly  upregulated  in  the  re-fed 
versus  the  -P  culture  at  48  h  (P<  0.001)  (Fig.  3). 

The  expression  of  nifH  was  examined  in  the  re-feed 
experiment.  In  the  -P  culture,  nifH  was  significantly  down¬ 
regulated  by  24  h  (P<  0.001)  relative  to  the  culture  re-fed 
with  P  (Fig.  4).  The  nifH  gene  continued  to  be  expressed  at 
the  same  level  in  the  cultures  where  P  was  added  (Fig.  4). 

Discussion 

Trichodesmium  populations  are  commonly  present  in  low- 
DIP  environments.  In  these  systems,  where  the  DOP  con¬ 
centration  typically  exceeds  that  of  DIP  in  the  euphotic 
zone  (Orrett  and  Karl,  1987;  Karl,  1999),  the  ability  to 
utilize  DOP  would  confer  a  distinct  growth  advantage  to 
marine  microbes.  Previous  studies  have  shown  that  cul¬ 
tures  and  field  populations  of  Trichodesmium  can  hydrol¬ 
yse  a  model  P  ester  (Stihl  eta!.,  2001;  Dyhrman  etai., 
2002).  Additionally,  it  has  been  suggested  that  Trichodes¬ 
mium  is  not  able  to  compete  with  co-occurring  plankton 
for  DIP,  and  that  it  uses  DOP  to  meet  its  P  demand  (Sohm 
and  Capone,  2006).  Linking  these  physiological  observa¬ 
tions  to  the  genes  that  control  these  activities  remains  a 
central  challenge  in  microbial  oceanography. 


Table  2.  Nucleotide  sequence  comparison  of  key  phosphorus-scavenging  genes  among  different  species  of  Trichodesmium. 


T.  tenue: 

T.  erythraeum 

T.  spiralis: 

T.  erythraeum 

T.  thiebautii: 

T.  erythraeum 

T.  tenue: 

T.  spiralis 

T.  tenue: 

T.  thiebautii 

T.  spiralis: 

T.  thiebautii 

psfS(1077)= 

99^ 

88 

88 

87 

88 

98 

sphX  (1086) 

98 

98 

99 

too 

99 

99 

phoA  (3528) 

98 

98 

98 

97 

96 

98 

phoX  (2358) 

99 

99 

97 

too 

98 

98 

a.  Size  of  putative  gene  in  T.  erythraeum  (bp). 

b.  %  iD  to  gene  fragment  sequenced.  The  psfS  fragment  is  -530  bp  for  all  species.  The  spbXfragment  is  -  600  bp  in  T.  tenue  and  T.  spiralis,  and 
is  400  bp  in  T.  thiebautii.  The  phoA  fragment  is  -1 .8  kb  in  T.  tenue,  240  bp  in  T.  thiebautii  and  202  bp  in  T.  spiralis.  The  phoX  fragment  is  1 36  bp 
in  aii  species. 
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Fig.  2.  Plots  of  APA  for  both  the  transfer  and  the  re-feed  experiments  (A  and  B),  and  expression  profiies  for  phoX  and  phoA  for  the  transfer 
experiment  (C  and  E),  and  the  re-feed  experiment  (D  and  F).  In  (A  and  B)  error  bars  indicate  standard  deviation  of  tripiicate  biological 
replicates.  In  (C-F)  error  bars  indicate  the  standard  error  of  the  fold  change  of  triplicate  measurements  between  the  sample  expression  and 
the  calibrator  expression  (i.e.  the  zero  time  point).  An  *  indicates  a  P-value  of  less  than  0.05  for  the  -P  versus  the  -i-P  or  re-feed  culture  at  a 
single  time  point  for  a  paired  f-test  (A  and  B)  or  a  pair-wise  fixed  reallocation  randomization  analysis  (C-F). 


As  highlighted  previously,  P  esters  dominate  the  high- 
molecular-weight  DOP  that  is  amenable  to  characteriza¬ 
tion  (Clark  eta!.,  1998),  and  P  esters  likely  represent  a 
significant  source  of  P  to  Trichodesmium  (Mulholland 
etal.,  2002;  Sohm  and  Capone,  2006).  Unfortunately,  the 
alkaline  phosphatases  that  hydrolyse  P  esters  can  be 
difficult  to  identify  by  sequence  homology  alone,  as  the 
primary  sequence  tends  to  be  very  heterogeneous 
(Moore  etal.,  2005).  Only  one  cyanobacterial  alkaline 
phosphatase  has  been  characterized:  the  atypical  alka¬ 
line  phosphatase  (phoA)  of  Synechococcus  PCC7942 
(Ray  etal.,  1991),  and  this  phoA  gene  has  very  little 
sequence  homology  to  that  of  other  characterized  phoA 
genes.  The  T.  erythraeum  genome  has  a  gene  similar  to 
this  atypical  alkaline  phosphatase  and  it  is  referred  to 


herein  as  phoA.  Another  alkaline  phosphatase,  phoX,  has 
recently  been  characterized  in  Pseudomonas  and  Pas- 
teurella  (Mends  etal.,  2006;  Wu  etal.,  2007),  and  it  con¬ 
sists  entirely  of  one  conserved  bacterial  domain  (COG 
3211).  Trichodesmium  has  two  genes  that  share  this  con¬ 
served  domain  (phoX and  phoX2).  Several  marine  cyano¬ 
bacteria  appear  to  have  this  gene,  including  Cyanothece 
ATCC  51142,  Microcystis  aeruginosa,  Synechococcus 
WHS  102,  Synechococcus  WH7803  and  Lyngbya  PCC 
8106,  which  has  two  copies  like  T.  erythraeum  IMS101 
(Sebastian  and  Ammerman,  2009).  Interestingly,  there  are 
no  apparent  phoX  genes  in  the  available  sequenced 
genomes  for  Prochlorococcus. 

The  putative  phoA  and  phoX  gene  sequences  were 
compared  between  all  cultured  Trichodesmium  species. 
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Fig.  3.  Expression  profiles  for  sphX  and  pstS  in  the  transfer  (A  and  C)  and  re-feed  (B  and  D)  experiments.  Error  bars  indicate  the  standard 
error  of  the  fold  change  of  triplicate  measurements  between  the  sample  expression  and  the  calibrator  expression  (i.e.  the  zero  time  point). 

An  *  indicates  a  P-value  of  less  than  0.05  for  the  -P  versus  the  -rP  or  re-feed  culture  at  a  single  time  point  using  a  pair-wise  fixed  reallocation 
randomization  analysis  (A-D). 


and  were  found  to  be  present  in  all  species  tested.  Over 
the  fragments  sequenced,  both  of  these  genes  were 
highly  conserved  within  the  genus,  suggesting  that  these 
genes  were  both  present  in  the  Trichodesmium  genome 
before  the  radiation  of  the  species. 


Fig.  4.  Expression  profiles  for  nifH  in  the  re-feed  experiment.  Error 
bars  indicate  the  standard  error  of  the  fold  change  of  triplicate 
measurements  between  the  sample  expression  and  the  calibrator 
expression  (i.e.  the  zero  time  point).  An  *  indicates  a  P-value  of 
less  than  0.05  for  the  -P  versus  re-feed  culture  at  a  single  time 
point  using  a  pair-wise  fixed  reallocation  randomization  analysis. 


The  phoA  and  phoX  genes  are  both  upregulated  by  P 
starvation,  but  phoX  showed  a  much  stronger  response 
than  phoA.  APAalso  showed  a  similar  pattern  of  induction 
to  the  gene  expression  pattern.  In  the  re-feed  experiment, 
phoA  and  phoX  expression  quickly  dropped  (by  4  h)  after 
P  was  added  back  to  a  P-starved  culture;  however,  APA 
decreased  more  slowly.  This  is  contrary  to  a  previous 
study  of  Trichodesmium  APA,  in  which  APA  decreased  by 
50%  or  more  within  48  h  after  P  was  added  back  to  a 
P-starved  culture  (Stihl  etal.,  2001).  Our  data  suggest 
that  the  alkaline  phosphatase  enzyme  may  persist  some¬ 
what  after  expression  of  the  gene  has  decreased,  and 
therefore  the  enzyme  is  not  actively  being  degraded  in  the 
time-course  of  our  study.  Although,  under  these  experi¬ 
mental  conditions,  we  cannot  rule  out  the  possibility  that 
additional  phosphatases  or  co-cultured  heterotrophs  are 
also  contributing  to  observed  activities.  In  this  context, 
measuring  phoA  or  phoX  gene  expression  may  provide  a 
more  instantaneous  picture  of  physiology,  whereas  APA 
may  reflect  a  longer  window  into  the  history  of  P  starva¬ 
tion.  This  is  an  important  consideration  given  the  common 
use  of  APA  as  a  physiological  indicator  of  P  limitation  in 
field  populations. 

This  work  represents  the  first  evidence  for  a 
P-regulated  marine  cyanobacterial  phoX.  The  difference 
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in  the  magnitude  of  expression  between  phoX  and  phoA 
in  Trichodemium  suggests  that  phoX  may  have  a  more 
significant  contribution  to  DOP  hydrolysis  than  phoA, 
although  there  need  not  be  a  linear  relationship  between 
expression  of  the  gene  and  activity  of  the  enzyme. 
Notably,  phoX  was  also  present  and  more  highly 
expressed  than  phoA  under  P  starvation  in  a  recent  study 
of  a  freshwater  thermophilic  Synechococcus  sp.  (Adams 
etai,  2008).  Given  the  presence  and  expression  ot  phoX 
in  several  cyanobacteria,  P  regulation  studies  in  other 
groups  would  assess  the  potential  Importance  of  this 
enzyme. 

If  the  phoX  and  phoA  alkaline  phosphatases  have  dif¬ 
ferent  substrate  specificity  and  cofactor  requirements,  it 
would  Increase  the  range  of  potential  DOP  compounds 
Trichodesmium  can  hydrolyse  and  the  range  of  trace 
metal  limitation  Trichodesmium  can  withstand.  The  char¬ 
acterization  work  done  with  PhoX  to  date  Indicates  a 
broad  substrate  specificity  including  some  diesters 
(Monds  etai.,  2006;  Wu  efa/.,  2007).  Additionally,  the 
PhoX  of  Pasteurella  does  not  require  Zinc  (Zn),  whereas 
proteins  of  the  PhoA  family  typically  require  Zn  as  a  cofac¬ 
tor  (Kim  and  Wyckoff,  1991;  Wu  etai.,  2007).  There  are 
several  regions  of  the  ocean,  such  as  the  Sargasso  Sea, 
where  both  P  concentrations  and  Zn  concentrations  are 
extremely  low  (Jakuba  etai.,  2008),  and  It  has  been 
hypothesized  that  the  Sargasso  Sea  Is  co-limited  by  P  and 
Zn  because  of  the  Zn  requirement  by  phosphatases 
(Shaked  etai.,  2006;  Jakuba  etai.,  2008).  Having  an 
alternative  alkaline  phosphatase  that  does  not  require  Zn 
as  a  cofactor  could  provide  a  growth  advantage  in  such 
areas.  Further  studies  into  the  substrate  specificity  and 
metal  requirements  of  the  PhoX  and  PhoA  in  Trichodes¬ 
mium  are  needed  to  address  this  observation. 

Previous  studies  of  the  presence  and  regulation  of  high- 
affinity  phosphate  binding  proteins  in  cyanobacteria  have 
provided  useful  Insights  Into  how  these  organisms  are 
adapted  to  low-P  environments  and  how  they  respond  to 
P  starvation  (Scanlan  etai.,  1997;  Moore  etai.,  2005; 
Martiny  etai.,  2006).  In  the  T.  erythraeum  genome  there 
appear  to  be  two  copies  of  genes  encoding  a  hIgh-affInIty 
P  binding  protein,  pstS  and  sphX.  Several  other  cyano¬ 
bacteria  have  multiple  copies  of  high-affinity  P  binding 
proteins,  for  example  In  freshwater  Synechocystis,  Gloeo- 
bacter,  Anabaena  PCC  7120,  and  In  marine  Synechococ¬ 
cus  WH81 02,  Synechococcus  WH7803,  Prochlorococcus 
MIT9313  as  well  as  the  marine  diazotroph  Crocosphaera 
WH8501  (Moore  etai.,  2005;  Dyhrman  and  Haley,  2006). 
Although  multiple  copies  are  common,  the  association  of 
pstS  and  sphX  with  pstCAB  ('responsible  for  P  transport 
into  the  cell)  and  the  gene  topologies  of  the  pstSCAB 
cluster  can  vary  substantially  between  genera  (Moore 
etai.,  2005;  Martiny  etai.,  2006).  Among  the  cyanobac- 
terial  genomes  currently  available,  the  topology  of  the 


Trichodesmium  Pst  system  Is  most  similar  to  that  of  Syn¬ 
echocystis,  although  the  sphX  and  pstS  genes  of  Tri¬ 
chodesmium  are  both  more  similar  to  the  pstS  gene  of 
Synechocytis  than  to  the  sphX  gene  of  Synechocystis. 

Given  the  marked  ecotype  variation  of  nutrient¬ 
scavenging  capacities  in  the  marine  cyanobacteria 
Prochiorococcus  and  Synechococcus  (Moore  et  a!.,  2005; 
Martiny  etai.,  2006),  we  were  interested  in  comparing  the 
sequence  and  presence  of  sphX  and  psfS  among  species 
of  Trichodesmium  (Moore  etai.,  2005).  All  species  tested 
contain  both  sphX and  pstS,  but  the  phylogenetic  relation¬ 
ship  Inferred  from  their  sequences  is  different.  There 
is  more  sequence  heterogeneity  among  the  species  in 
pstS  than  sphX.  The  sequence  identity  for  pstS  divides 
the  species  into  two  distinct  clusters:  T.  erythraeum  and 
T.  tenue  In  one  cluster  and  T.  thiebautii  and  T.  spiratis  In 
another.  This  differs  from  the  phylogeny  Inferred  using  the 
Internal  transcribed  spacer  sequence,  which  places 
T.  erythraeum  In  one  clade  and  T.  tenue,  T.  spiralis,  and 
T.  thiebautii  \n  another  (Orcutt  etai.,  2002),  and  the  rela¬ 
tionship  Inferred  from  the  sequence  of  a  functional  gene 
sequence,  hetR,  which  places  T.  spiralis  and  T.  thiebautii 
In  one  cluster  and  T.  tenue  and  T.  erythraeum  each  In  their 
own  clusters  (Lundgren  etai,  2005).  Given  the  sequence 
heterogeneity  in  pstS  found  here,  further  analysis  of  the 
pstSCAB  topology  and  expression  patterns  Is  warranted 
In  the  other  cluster. 

The  sphX  gene  was  regulated  by  P  supply  In  the  trans¬ 
fer  experiment,  while  pstS  was  not.  The  induction  time  for 
sphX  was  longer  than  what  was  observed  for  phoA  and 
phoX',  however,  It  Is  difficult  to  resolve  significant  expres¬ 
sion  changes  for  sphX  relative  to  these  genes  because  of 
the  lower  expression  levels.  This  Is  a  relatively  long  lag 
before  Induction  of  sphX  in  comparison  with  other  studies 
of  cyanobacterial  pstS/sphX  expression  in  the  unicellular 
Prochiorococcus  and  Synechocystis,  which  tend  to 
respond  within  48  h  (Suzuki  etai.,  2004;  Martiny  etai., 
2006).  However,  it  has  been  observed  that  P  can  be 
adsorbed  to  the  outside  of  Trichodesmium  cell  walls 
(Sahudo-Wilhelmy  etai.,  2004),  and  this  P  may  not  be 
removed  by  the  washing  procedure  used  in  this  experi¬ 
ment.  This  could  delay  the  onset  of  P  starvation,  as  would 
Internal  stores  of  P,  such  as  polyphosphate.  Additionally, 
Trichodesmium  Is  relatively  slow-growing  and  may  not 
sense  and  respond  to  Its  environment  as  quickly  as  other 
cyanobacteria. 

In  the  re-feed  experiment,  sphX  was  downregulated  in 
response  to  the  P  addition,  but  pstS  was  not.  In  fact,  pstS 
was  upregulated  at  the  48  h  time  point.  The  global  down- 
regulation  of  genes  In  response  to  starvation,  and  recov¬ 
ery  from  P  starvation  in  the  re-fed  culture  may  be  driving 
this  change  in  psfS  expression.  While  there  was  no  clear 
P  regulation  of  pstS  In  these  experiments,  it  is  possible 
that  a  slightly  different  time-course,  or  more  sensitive 
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tools,  would  reveal  subtle  regulation  patterns  that  we  were 
unable  to  detect. 

Differential  regulation  of  genes  encoding  phosphate 
binding  proteins  has  been  observed  In  other  cyanobacte¬ 
ria  (Moore  etal.,  1998;  Suzuki  etal,  2004),  but  never 
between  two  such  genes  that  appear  to  be  In  a  single 
operon.  For  example,  Synechocystis  has  two  non¬ 
identical  copies  of  the  entire  Pst  system  (Psti  and  Pst2). 
Both  copies  are  regulated  by  P,  but  one,  Psti ,  acts  as  a 
rapid  response  to  P  starvation  while  Pst2  is  upreguiated 
later  under  P-starvation  conditions  (Suzuki  etal.,  2004). 
Prochlorococcus  MIT9313  also  has  two  copies  of  pstS, 
one  regulated  by  P  and  one  not,  but  neither  are  associ¬ 
ated  with  pstCAB  as  is  the  case  in  T.  erythraeum  (Martiny 
etal.,  2006).  Given  the  location  of  the  pstS  gene  relative 
to  the  rest  of  the  Pst  system,  it  is  possible  that  pstCAB  are 
co-transcribed  with  pstS  and  this  suggests  that  pstCAB 
are  expressed  but  not  regulated  by  P.  With  the  apparent 
lack  of  the  constitutively  expressed  Pit  system,  it  is  pos¬ 
sible  that  the  PstSCAB  system  is  fulfilling  that  role  for  P 
uptake  during  replete  growth  conditions,  with  sphX 
upreguiated  under  low  P.  This  is  consistent  with  the  posi¬ 
tion  of  the  Pho  box  (see  below)  and  with  studies  on 
Trichodesmium  P  uptake  kinetics,  where  under  P  limita¬ 
tion  the  Vmax  for  P  uptake  increases,  but  the  Km  remains 
constant  relative  to  replete  cultures  (Fu  etal.,  2005). 

To  evaluate  the  diagnostic  effectiveness  of  the  pres¬ 
ence  of  a  Pho  box  for  identifying  P-regulated  genes,  we 
summarized  ail  the  cyanobacteria  that  have  been 
examined  for  Pho  boxes  (Su  etal.,  2007)  and  for  which 
there  are  P  regulation  data  (Table  1).  In  most  species, 
P-regulated  genes  are  preceded  by  a  Pho  box,  except  in 
Prochlorococcus  marinus  MiT9313  for  which  Su  and  col¬ 
leagues  (2007)  did  not  identify  a  Pho  box  preceding  the 
P-regulated  pstS  or  phoA  (Martiny  etal,  2006;  Su  etal., 
2007).  Also,  while  there  was  no  Pho  box  preceding  phoA 
in  T.  erythraeum,  this  gene  was  still  regulated  by  P.  The 
presence  of  a  sequence  resembling  a  Pho  box  preceding 
this  gene  suggests  that  there  is  possibly  more  heteroge¬ 
neity  in  this  binding  sequence  in  Trichodesmium  than 
predicted,  in  Synechocystis,  sphX  and  pstSCAB  are 
separated  by  a  small  gap,  but  both  sphX  and  pstS  are 
preceded  by  a  Pho  box  (Su  etal.,  2007).  In  T.  erythraeum, 
sphX  and  pstSCAB  are  also  separated  by  a  small  gap, 
and  only  sphX  is  preceded  by  a  putative  Pho  box  (Su 
etal.,  2007).  The  absence  of  a  second  Pho  box  is  con¬ 
sistent  with  the  lack  of  P  regulation  in  pstS.  Therefore, 
examining  a  genome  for  Pho  box  sequences  is  a  useful 
first  tool  for  identifying  potential  genes  involved  in  P  acqui¬ 
sition,  but  the  regulation  and  activity  of  these  genes  still 
needs  to  be  confirmed. 

To  examine  the  impact  of  P  starvation  on  N2  fixation, 
nifH  (nitrogenase  reductase)  expression  was  assayed  as 
a  proxy  for  N2  fixation  rate  (Chen  et  a!.,  1 998;  El-Shahawy 


eta!.,  2003).  In  the  re-feed  experiment,  nifH  expression 
was  significantly  downregulated  in  the  P-starved  cultures, 
and  therefore  P  starvation  may  lead  to  a  decrease  in  nifH 
expression  and  a  reduction  in  N2  fixation  rate.  This  has 
important  implications  for  the  field,  and  suggests  that 
increased  expression  of  sphX,  phoX or  phoA,  indicates  a 
decrease  in  N2  fixation  due  to  P  starvation  in  situ.  When 
re-fed  with  P,  nifH  expression  remained  constant  relative 
to  the  start  of  the  experiment  while  the  P-regulated 
genes  were  repressed,  suggesting  that  the  link  between 
decreased  N2  fixation  rates  and  P  starvation  is  not  as  tight 
as  has  been  seen  with  iron  starvation  (Shi  etal.,  2007; 
Kupper  etal.,  2008).  A  more  detailed  time-course  study 
with  N2  fixation  assays  would  further  resolve  molecular- 
level  linkages  between  P  starvation  and  changes  in  N2 
fixation  rate. 

Here  we  studied  the  interplay  among  P  starvation,  APA 
and  gene  expression  patterns.  The  pstS,  sphX,  phoA  and 
phoX  genes  are  present  in  all  examined  species  of  Tri¬ 
chodesmium.  Additionally,  the  expression  of  phoA,  phoX 
and  sphX  were  all  responsive  to  P  supply,  while  pstS  was 
not.  These  data  suggest  that  the  expression  of  sphX, 
phoA  and  phoX could  provide  a  useful  marker  for  P  physi¬ 
ology  in  natural  populations.  Combining  these  measure¬ 
ments  with  APA  would  also  provide  insights  into  the 
nutritional  history  of  this  organism.  Together  with  mea¬ 
surements  of  n/TH  expression  or  N2  fixation  rates,  we  can 
identify  the  importance  and  prevalence  of  P  starvation, 
and  decreases  in  N2  fixation  due  to  P  starvation  in  the 
environment.  Additionally,  we  demonstrate  for  the  first 
time  the  P  regulation  of  a  putative  phoX  in  a  marine 
cyanobacterium.  The  presence  of  a  P-regulated  phoX  in 
Trichodesmium  could  provide  a  means  to  access  the  P 
ester  pool  without  an  increased  Zn  requirement.  In  olig- 
otrophic  regimes,  where  P  esters  are  the  dominate  P  pool 
and  Zn  concentrations  are  low  (Jakuba  etal.,  2008), 
having  a  PhoX-type  alkaline  phosphatase  could  allow 
Trichodesmium  to  occupy  a  unique  niche  relative 
to  cyanobacteria  solely  using  a  PhoA-type  alkaline 
phosphatase. 

Experimental  procedures 

Culture  conditions 

Trichodesmium  cultures  were  grown  in  BMP  medium  with 
Sargasso  seawater  as  previously  described  (Webb  etal., 
2001 ).  Species  examined  include  Trichodesmium  erythraeum 
IMS101,  T.  tenue,  T.  thiebautii  and  T.  spiralis,  obtained 
from  the  culture  collection  of  John  Waterbury.  All  cultures 
were  unialgal  but  not  axenic.  Cultures  were  grown  at  25 °C 
on  a  shaker  table  with  daily  cycles  that  consisted  of 
1h  at  13.5  |iE  s ',  10  h  at  33.7  |iE  s ',  1h  at 
1 3.5  pE  m^^  s  ’  and  1 2  h  of  dark.  To  examine  the  regulation  of 
P  acquisition  genes  in  the  transfer  experiment,  1  I  of  expo¬ 
nentially  growing  cells  was  divided  in  half  and  gently  filtered. 
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Table  3.  Primer  sequences  and  product  sizes  for  PCR  and  qPCR  reactions. 


Gene 

Sequencing  primers  5'-3' 

Annealing 

Temperature  (°C) 

Product 

size 

Efficiency 

phoA 

ATGCGTGGGGACTTAACAGTAA 

TCTAATCACAAAATCATCTGTTGTGAGAG 

61.2 

3.5  kb 

N.A. 

pstS 

ACAAGCACAAACTAAAACCAG 

GACGAATCAGCAGTGACAAG 

58.5 

1.5  kb 

N.A. 

sphX 

AAATTAGTATCGTTGCCTAAAT 
TTTGTGCTTATACATTTATTATC 
qPCR  primers  5'-3' 

59 

1.5  kb 

N.A. 

phoA 

AGATGATGATGGAGAAAATGAAAG 

CGTCTTCCCTGGTGAGGAG 

60.7 

131  bp 

105% 

phoX 

TGGTGAAAAGTATGGCTGGATGG 

GCTGCTAACGGTTTACCTGCTTC 

60.7 

135  bp 

99% 

pstS 

TGGAGTTTTCACGAAGCATTTGAG 

GCAGTTACACCATCATTCTTGGC 

60.7 

120  bp 

98% 

sphX 

TCAGTTGGTAGTGGTGCTGG 

CGAGAGACCTTGGCTATTTCAC 

60.7 

104  bp 

103% 

nifH 

AGTGTGTTGAGTCTGGTGGTCCTG 

TGTTTTCACGGATAGGCATTGC 

60.7 

174  bp 

101% 

rotA 

CCAAGGAGGTTGCCCAAAAGG 

GAACGCCATCTAAATGAGGTTGAG 

60.7 

175  bp 

97% 

Annealing  temperatures  for  PCR  reactions  and  primer  efficiencies  are  as  indicated.  If  no  sequencing  primer  is  listed,  then  qPCR  primers  were  used 
for  sequencing. 

N.A.,  not  applicable. 


washed  three  times  in  media  with  or  without  8  |iM  H3PO4 
(+P  or  -P),  then  re-suspended  in  500  mi  of  -i-P  or  -P  media 
respectiveiy.  The  resuiting  ceii  suspensions  were  spiit  into 
repiicate  500  mi  poiycarbonate  fiasks  with  100  mi  of  cuiture  in 
each.  Whoie  fiasks  were  harvested  at  0,  24,  72,  168  and 
216  h.  For  each  time  point,  -i-P  and  -P  fiasks  were  harvested 
one  at  a  time  starting  approximateiy  2  h  after  the  mid-point  of 
the  iight  cycie.  For  this  experiment,  we  couid  not  assay  nifH 
expression  due  to  the  diei  periodicity  of  the  expression  of  this 
gene  and  differences  in  harvest  time  between  sampies.  To 
examine  the  downreguiation  of  P  acquisition  genes  in  the 
re-feed  experiment,  ceiis  were  grown  without  P  for  10  days, 
and  P  was  added  back  to  two  of  the  fiasks  to  a  concentration 
of  8  pM.  Fiasks  were  subsampied  at  0,  4,  24  and  48  h.  For 
each  time  point,  fiasks  were  harvested  in  paraiiei  exactiy  2  h 
after  the  mid-point  of  the  iight  cycie,  aiiowing  for  n/ff-/ expres¬ 
sion  assays. 

Genome  annotation 

P-scavenging  genes  were  identified  using  the  integrated 
microbiai  genomes  portai  at  the  Joint  Genome  institute-US 
Department  of  Energy  website  (http://img.jgi.doe.gov/) 
and  the  Nationai  Center  for  Biotechnoiogy  information 
GenBank  database  (http://www.ncbi.nim.nih.gov/),  based  on 
homoiogy. 

DNA  amplification  and  sequencing 

Ceiis  were  coiiected  by  fiitration  onto  a  5  pm  poiycarbonate 
fiiter  and  stored  at  -20  °C  untii  extraction.  DNA  was  extracted 
using  the  instagene  Matrix  (Bio-Rad,  Flercuies,  CA)  accord¬ 
ing  to  the  manufacturer’s  instructions  or  as  described  eise- 
where  (Ehrenreich  eta!.,  2005).  PCR  primers  were  designed 
to  either  ampiify  fragments  or  the  entire  putative  sphX,  pstS, 


phoX  and  phoA  genes.  Anneaiing  temperatures  were  empiri- 
caiiy  determined  for  each  primer  pair  (Tabie  3).  Each  PCR 
reaction  consisted  of  a  5  min  denaturation  step  at  95 °C, 
foiiowed  by  35  cycies  of  1  min  at  95 °C,  1  min  at  the  anneaiing 
temperature,  1  min  per  500  bp  at  72 °C  and  a  finai  extension 
of  10  min  at  72‘G,  in  a  Bio-Rad  iCycier  (Flercuies,  CA).  PCR 
amplification  was  done  using  0.5  pi  PfuTurbo  DNA  poly¬ 
merase  (Stratagene,  La  Jolla,  Ca),  2.5  pi  of  DNA  template, 
2.5  pi  of  2  mM  dNTP,  2.5  pi  of  1 0x  Pfu  buffer,  50  pmol  of  each 
primer  and  sterile  water  to  a  final  volume  of  25  pi.  PCR 
products  were  gel-extracted  using  the  QIAquick  gel  extrac¬ 
tion  kit  (Qiagen,  Valencia,  CA),  and  direct-sequenced  at  the 
Josephine  Bay  Paul  Center  of  the  Marine  Biological  Labora¬ 
tory  (Woods  Hole,  MA),  MWG  Biotech  (Hunstville,  AL)  or  at 
the  University  of  Maine  (Orono,  ME)  according  to  the  facility’s 
protocols.  Using  the  genome  of  T.  erythraeum  IMS101  as  a 
guide,  internal  sequencing  primers  were  designed  for  any 
PCR  product  sequenced  larger  than  500  bp  in  order  to  cover 
the  entire  product  at  least  once  on  each  strand.  Sequences 
were  edited  using  Sequencher  (Gene  Codes,  Ann  Arbor,  Ml) 
and  verified  manually.  Alignments  were  done  in  MacVector 
(Cary,  NC).  The  sequences  have  been  deposited  in  General 
Bank  (Accession#  FJ602760-602771). 


RNA  isolation  and  cDNA  synthesis 

Cells  were  collected  via  filtration  on  5  pm  polycarbonate 
filters  and  immediately  stored  in  liquid  nitrogen  until  extrac¬ 
tion.  RNA  was  extracted  using  RNeasy  kit  (Qiagen,  Valencia, 
CA)  according  to  the  manufacturer’s  instructions  using  the 
bead  beating  protocol.  The  RNA  was  then  treated  with  Turbo 
DNase  (Ambion,  Austin,  TX)  following  the  rigorous  protocol  to 
remove  genomic  DNA.  RNA  was  quantified  spectrophometri- 
cally,  and  a  total  of  60  ng  of  RNA  was  transcribed  into  cDNA 
using  the  iScript  cDNA  Synthesis  kit  (Bio-Rad,  Hercules,  CA). 
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A  second  reaction  was  performed  for  each  sample  containing 
no  reverse  transcriptase  to  serve  as  controls  for  genomic 
contamination  in  subsequent  analyses. 

Quantitative  PCR 

Primers  were  designed  using  MacVector  (Cary,  NC).  For 
each  biological  replicate,  triplicate  reactions  were  run  using 
the  iQ  SYBR  Green  SuperMix  (Bio-Rad,  Hercules,  CA)  as 
follows:  12.5  pi  SYBR  Green  SuperMix,  200  nM  forward  and 
reverse  primers,  2  pi  of  diluted  cDNA  (the  rough  equivalent  of 
1.2  ng  of  total  RNA),  and  sterile  water  to  a  final  volume  of 
25  pi.  The  Quantitative  PCR  (qPCR)  programme  consisted  of 
95 °C  for  5  min,  followed  by  45  cycles  of  95 °C  for  30  s,  60.7°C 
for  30  s,  then  one  cycle  of  95 °C  for  1  min  in  an  iCycler  iQ 
(Bio-Rad,  Hercules,  CA).  Melt  curves  for  each  reaction  were 
performed  to  test  the  quality  of  the  product.  For  each  sample 
and  primer  pair  no  reverse  transcriptase  controls  and  no 
template  controls  were  run  to  test  for  DMA  contamination  and 
specificity  of  amplification.  These  negative  controls  did  not 
typical  amplify,  but  in  the  rare  cases  where  there  was  any 
amplification  the  measured  template  concentration  was  at 
least  two  orders  of  magnitude  less  than  the  sample  template 
concentration. 

The  comparative,  or  AACt,  method  was  used  to  normalize 
the  expression  of  each  gene.  The  amplification  efficiency  and 
validation  of  the  AACt  method  were  tested  for  each  primer  set 
on  cDNA.  Amplification  efficiency  was  determined  to  be 
acceptable  if  it  fell  between  90%  and  110%  with  an  ^  of 
0.98  or  higher  (Table  3).  Each  primer  pair  was  considered 
acceptable  if  the  ACt  of  the  two  slopes  was  less  than  0.1. 
Samples  were  endogenously  normalized  to  the  constitu- 
tively  expressed  rotA  Cencoding  a  peptidyl-prolyl  cis-trans 
isomerase)  (Hassidim  etal.,  1992;  Thellin  etai.,  1999).  The 
Ct  of  the  reference  gene  did  not  vary  with  treatment,  further 
validating  its  application  as  a  reference  gene.  Fold  changes 
were  calculated  using  the  Relative  Expression  Software 
Tool,  which  accounts  for  differences  in  efficiency  in  its 
fold  change  calculation  (http://www.gene-quantification.de/ 
download.html).  For  each  experiment  a  single  To  time  point 
sample  was  used  as  a  calibrator  for  calculating  fold  changes. 
The  P-values  were  calculated  using  a  pair-wise  fixed  reallo¬ 
cation  randomization  analysis  comparing  the  -i-P  to  the  -P 
samples  at  each  time  point  (PfaffI  etal.,  2002). 

Chlorophyll  a 

Between  5  and  25  ml  of  culture  were  collected  onto  GF/F 
filters  and  frozen  at  -20 °C.  Filters  were  extracted  in  90% 
acetone  and  incubated  overnight  at  -20 °C.  Chlorophyll  a 
fluorescence  was  measured  on  an  aquafluor  handheld  fluo- 
rometer  and  analysed  according  to  Arar  and  Collins  (1997). 

Bulk  APA 

For  APA  assays,  triplicate  samples  at  each  time  point  were 
taken.  These  consisted  of  5-20  ml  of  culture  filtered  onto  a 
5  pm  polycarbonate  filter  and  frozen  at  -20 °C.  Samples  were 
processed  as  described  elsewhere  (Dyhrman  and  Rutten- 
berg,  2006).  Briefly,  lOpM  6,8-difluoro-4-methylumbelliferyl 


phosphate  was  added  to  each  sample  in  a  Petri  dish  and 
fluorescence  was  measured  on  a  Fluostar  plate  reader  every 
5-20  min  for  five  time  points.  Previous  kinetics  experiments 
found  the  10  pM  substrate  concentration  to  be  saturating 
(data  not  shown).  Standard  curves  were  generated  for  each 
sample  run  using  6,8-difluoro-7-hydroxy-4-methylcoumarin. 
The  -i-P  and  -P  samples  at  each  time  point  were  determined 
to  be  statistically  significantly  different  using  a  paired  f-test. 
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Chapter  5 


Molecular  evidence  for  phosphorus  limitation  of  Trichodesmium  growth  and 
nitrogen  fixation  in  the  North  Atlantic 


Abstract 

The  nitrogen  (N)  fixing  cyanobacterium,  Trichodesmium  has  a  significant  impact 
on  the  N  and  carbon  (C)  cycles  of  the  tropical  and  subtropical  oceans.  Phosphorus  (P) 
concentrations  are  very  low  in  these  regions,  and  therefore  P  has  been  suggested  to 
constrain  N2  and  CO2  fixation  by  this  genus.  However,  there  is  little  direct  evidence  of  P 
limitation  of  Trichodesmium  or  its  affect  on  N2  fixation  in  these  tropical  oceans.  By 
combining  molecular  techniques  with  culture  experiments  we  can  probe  the  interplay 
between  nutrient  limitation,  gene  expression,  and  physiological  responses  such  as 
decreased  N2  fixation  rates  in  culture  and  in  the  field.  This  study  identifies  P  limitation 
levels  that  are  indicative  of  decreases  in  N2  fixation  rate  in  Trichodesmium  by  measuring 
N2  fixation  and  several  metrics  of  P  physiology  (alkaline  phosphatase  activity,  cellular 
ratios  of  C:P  and  phoX  gene  expression)  in  culture.  Expression  of  the  phoX  gene,  which 
encodes  an  alkaline  phosphatase,  was  calibrated  to  decreases  in  N2  fixation  in  semi- 
continuous  cultures  of  T.  erythraeum  IMS  101  grown  at  different  P-limited  growth  rates. 
The  expression  of  phoX,  measured  in  field  populations  of  Trichodesmium  from  the  North 
Atlantic  and  South  Pacific,  was  significantly  correlated  to  inorganic  phosphate 
concentration.  Additionally,  expression  levels  of  phoX  in  the  North  Atlantic  were 
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indicative  of  a  40-70%  reduction  in  N2  fixation  due  to  P  limitation.  These  data  indicate 


that  P  can  limit  growth  and  N2  fixation  in  Trichodesmium  from  the  North  Atlantic. 

Introduction 

Understanding  the  controls  on  primary  production  is  a  critical  area  of  research  in 
oceanography.  Marine  diazotrophs,  such  as  Trichodesmium,  are  important  contributors 
to  primary  production  because  they  fix  both  atmospheric  carbon  (C)  and  nitrogen  (N). 
Trichodesmium  often  dominates  in  tropical  and  subtropical  oceans  where  approximately 
50%  of  the  oceans  primary  production  occurs  (Capone  et  ah,  1997;  Karl,  2002).  In  these 
regions  the  input  of  N  from  diazotrophs  accounts  for  an  estimated  50%  of  the  “new” 
nitrogen  to  these  systems  (Karl,  2002).  Further,  it  is  this  “new”  nitrogen  that  supports 
export  production  in  oligotrophic  regions  (Capone,  2001).  Therefore,  understanding  the 
controls  on  growth  and  N2  fixation  of  Trichodesmium  is  critical  to  understanding  the 
constraints  on  primary  productivity  in  oligotrophic  oceans. 

Iron  (Fe)  and  phosphorus  (P)  are  generally  considered  the  two  most  important 
nutrients  in  limiting  N2  fixation.  Culture  work,  field  correlations,  and  qualitative 
molecular  assays  have  repeatedly  suggested  that  Trichodesmium  can  be  Fe-limited 
(Berman-Frank  et  ah,  2001;  Kustka  et  ah,  2003;  Moore  et  ah,  2009)  P-1  i  mi  ted  (Wu  et  ah, 
2000;  Sanudo-Wilhelmy  et  ah,  2001;  Dyhrman  et  ah,  2002;  Sohm  et  ah,  2008;  Hynes  et 
ah,  2009),  or  Fe  and  P  co-limited  (Mills  et  ah,  2004).  The  bulk  of  these  studies  suggest 
that  Trichodesmium  or  N2  fixation  is  more  often  Fe-limited  in  the  Pacific  and  South 
Atlantic,  and  P-limited  in  the  North  Atlantic.  These  conclusions  are  based  on  studies 
using  a  variety  of  metrics  to  determine  nutrient  physiology,  and  the  results  of  these 
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experiments  can  vary  greatly  depending  on  the  time  scale  considered  and  the  process 
being  measured. 

Traditional  metrics  for  determining  nutrient  limitation,  such  as  incubation 
experiments,  chemical  analysis,  and  mathematical  modeling  have  been  extremely  useful, 
but  they  do  not  allow  for  the  determination  of  the  nutritional  status  of  a  single  component 
in  a  natural  assemblage.  Enzyme-labeled  fluorescence  (ELF)  alkaline  phosphatase 
activity  (APA)  assays  do  allow  for  cell-specific  identification  of  P  physiology,  but  is 
qualitative  (Dyhrman  et  ah,  2002).  There  is  a  need  for  new  techniques  that  can  probe  the 
physiology  of  a  single  group  of  interest  quantitatively.  Quantitative  gene  expression  is 
just  such  an  assay.  Previous  work  has  shown  that  Trichodesmium  regulates  phoX,  which 
encodes  an  alkaline  phosphatase,  in  response  to  P  physiology  (Orchard  et  ah,  2009). 
Therefore,  the  expression  of  this  gene  in  the  field  can  be  a  genus-specific  method  to 
quantitatively  assess  P  physiology.  Additionally,  the  abundance  of  the  phoX  transcript 
decreases  rapidly  upon  recovery  from  P  limitation  in  culture  (within  hours),  unlike  the 
activity  of  the  enzyme  it  encodes  (APA)  which  can  take  days  to  decrease  (Orchard  et  ah, 
2009).  Because  of  its  rapid  response  and  turnover,  phoX  expression  is  a  more 
instantaneous  metric  of  P  physiology  than  APA. 

This  study  used  a  semi-continuous  culture  technique  to  link  P-limited  growth  rate 
to  changes  in  N2  fixation  rate,  APA,  cellular  C:P  ratios,  and  phoX  expression  in 
laboratory  cultures.  The  expression  of  phoX  was  then  measured  in  the  both  North  Atlantic 
and  South  Pacific  Trichodesmium  populations,  providing  a  direct  and  quantitative  assay 
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of  P  limitation  of  Trichodesmium  N2  fixation.  To  our  knowledge,  this  is  the  first 
quantitative  genus-specific  detection  of  P  limitation  in  Trichodesmium  field  populations. 

Methods 

Culture  experiments 

Trichodesmium  erythraeum  IMS  101  cultures  were  grown  in  250  mL  flasks  at 
25°C  on  a  shaker  table  with  daily  cycles  that  consisted  of  1  hr  at  13.5  pEin  m'^  s'*,  10  hrs 
at  33.7  pEin  m'  s'  ,  1  hr  at  13.5  pEin  m'  s'  ,  and  12  hrs  of  dark.  Cultures  were  grown  in 
duplicate  on  RMP  media  as  described  in  Webb  et  al.  (2001)  using  a  phosphate 
concentration  of  2  pM.  Once  cultures  reached  exponential  growth  they  were  diluted 
daily.  Two  separate  semi-continuous  experiments  were  preformed,  one  to  measure  N2 
fixation,  APA,  C:P  ratio,  which  was  diluted  at  0.1,  0.15  or  0.25  d'*.  The  second  measured 
phoX  and  nifH  expression,  and  the  cultures  were  diluted  at  0.1,  0.2  or  0.25  d'*.  Cultures 
were  considered  entrained  at  a  specific  growth  rate  when  in-vivo  chlorophyll  a 
fluorescence  was  stable  for  seven  days.  These  cultures  were  unialgal,  but  not  axenic. 
However,  biomass  was  dominated  by  Trichodesmium. 

Field  Samples 

Samples  were  taken  during  three  cruises  to  the  North  Atlantic  and  one  in  the 
South  Pacific  (Eigure  1;  Table  1).  Trichodesmium  colonies  were  collected  from  the  near 
surface  (roughly  within  the  top  20m)  using  a  hand  held  130  jim  net.  Colonies  were  split 
by  morphotype  and  only  colonies  with  a  raft  morphology  (parallel  trichomes)  were  used 
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in  further  analyses.  Single  colonies  were  picked  into  0.2  /^m-filtered  local  surface 
seawater,  and  then  transferred  into  fresh  0.2  /^m-filtered  water  to  reduce  contamination  of 
closely  associated  organisms.  Colonies  were  then  separated  for  further  analyses  as 
indicated  below.  To  the  extent  possible,  samples  were  collected  at  similar  times  of  the 
day  to  avoid  possible  diel  variability  in  gene  expression  patterns  (Table  1).  Water 
samples  were  collected  from  a  Niskin  bottle  on  a  conductivity,  temperature  and  depth 
(CTD)  rosette  from  a  depth  of  5  m,  and  were  processed  to  determine  nutrient 
concentrations  as  described  below. 

Enzyme  activity  assays 

For  APA  assays,  3-5  field  colonies  or  5-20  ml  of  culture  were  filtered  onto  a  5/^m 
polycarbonate  filter.  Samples  were  processed  as  described  elsewhere  (Dyhrman  and 
Ruttenberg,  2006).  Briefly,  10  /^mol  L'^  6,8-difluoro-4-methylumbelliferyl  phosphate 
(di-MUF-P)  (Invitrogen)  was  added  to  each  sample  in  a  petri  dish  with  artificial  seawater 
containing  no  added  P,  and  fluorescence  was  measured  on  a  Fluostar  Optima  plate  reader 
(BMG  Labtech)  every  5-20  minutes  for  5  time  points,  within  the  linear  range  of  the 
assay.  Previous  kinetics  experiments  found  the  10  //mol  L'^  substrate  concentration  to  be 
saturating  (data  not  shown).  Standard  curves  were  generated  for  each  assay  using  6,8- 
difluoro-7 -hydroxy-4-methylcoumarin  (di-MUF) . 

N2  fixation  was  measured  using  the  acetylene  reduction  assay  as  described  in 
Capone  and  Montoya  (2001).  Incubations  were  conducted  with  30  ml  of  culture  in  75  ml 
(total  volume)  Nalgene  bottles  with  septa  lids.  To  each  sample,  6  ml  of  acetylene  was 
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added  to  the  headspace  at  the  mid-point  of  the  light  cycle  (when  Trichodesmium  N2 
fixation  should  be  highest).  The  ethylene  concentration  in  the  headspace  was  measured 
in  duplicate  every  hour  using  a  GC-8A  gas  chromatograph  (Shimadzu  Norwell,  MA).  A 
ratio  of  4: 1  C2H2:N2  was  used  to  convert  to  N2  fixation  rates.  A  two-tailed  t  test  was  used 
to  determine  if  rates  were  statistically  different. 

Chlorophyll  a 

Between  5-25  ml  of  culture  was  collected  onto  GF/F  filters  and  frozen  at  -20°C. 
Filters  were  extracted  in  90%  acetone  and  incubated  overnight  at  -20°C.  Chlorophyll  a 
fluorescence  (chi  a)  was  measured  on  an  aquafluor  handheld  fluorometer  and  analyzed 
according  to  Arar  and  Collins  (1997). 

Nutrients 

Approximately  15  ml  of  culture  from  the  semi-continuous  experiment  was 
collected  on  pre-combusted  (500°C  for  5  hours)  25mm  Whatman  GF/F  and  dried  at  60°C 
for  particulate  C  and  P  measurements.  Particulate  C  was  measured  on  Flash  EAl  1 12 
CHN  Analyzer  using  a  Dynamic  Flash  Combustion  technique.  Intracellular  P  (Pic)  and 
total  cellular  P  (Ptc)  were  compared  by  washing  or  not  washing  with  oxalate  reagent  as 
in  (Sanudo-Wilhelmy  et  ah,  2004).  Samples  were  then  combusted  for  1  hour  at  500°C. 
They  were  then  incubated  at  104 °C  for  2  hours  with  0.4  ml  of  1  N  HCl  and  2  ml  of  DI 
water,  and  measured  according  the  method  of  Murphy  et  al.  (1962).  Samples  were  read 
on  a  UV-1700  Spectrophometer  (Shimadzu,  Norwell  MA)  at  885nm.  A  two-tailed  t  test 
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was  used  to  determine  if  ratios  were  statistically  different. 

Soluble  reactive  phosphorus  (SRP)  from  a  depth  of  5  m  was  measured  according 
to  the  MAGIC-SRP  method  (Karl  and  Tien,  1992),  for  samples  from  the  Atlantic  or  using 
the  molybdenum  blue  procedure  (Bernhardt  and  Wilhelms,  1967)  for  samples  from  the 
Pacific.  For  station  1 1  on  the  SPAC  cruise,  there  was  no  SRP  sample  taken,  but  it  is 
estimated  as  54  nM  as  the  average  of  the  SRP  concentration  from  the  two  closest  stations 
(Table  1).  SRP  concentrations  measured  here  are  referred  to  as  the  dissolved  inorganic 
phosphorus  (DIP)  concentration. 

RNA  Isolation  and  cDNA  Synthesis 

Field  colonies  or  cultures  were  collected  via  filtration  on  5  pm  polycarbonate 
filters,  frozen  in  liquid  nitrogen  as  quickly  as  possible  (typically  around  2  minutes  for 
cultures  and  15-30  minutes  for  field  samples),  and  stored  in  liquid  nitrogen  until  they 
were  extracted.  RNA  was  extracted  using  an  RNeasy  kit  (Qiagen,  Valencia,  CA) 
according  to  the  manufacturer’s  instructions  using  the  bead  beating  protocol.  The  RNA 
was  then  treated  with  Turbo  DNase  (Ambion,  Austin,  TX)  following  the  rigorous 
protocol  to  remove  genomic  DNA.  RNA  from  field  samples  were  precipitated  with 
100%  ethanol  and  stored  at  -80°C  overnight.  Field  RNA  samples  were  then  spun  for  30 
minutes  at  14000  rpm  and  the  ethanol  was  removed  form  the  RNA  pelleted,  and  the  pellet 
was  re-suspended  in  RNase  free  water.  All  RNA  was  quantified  spectrophometrically, 
and  a  total  of  60  ng  of  RNA  for  culture  samples  and  120  ng  for  field  samples  was 
transcribed  into  cDNA  using  the  iScript  cDNA  Synthesis  kit  (Bio-Rad,  Hercules,  CA).  A 
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second  reaction  was  performed  for  each  sample  containing  no  reverse  transcriptase  to 
serve  as  controls  for  genomic  contamination  in  subsequent  analyses. 

Quantitative  PCR  (qPCR) 

Quantitative  PCR  was  performed  following  the  protocol  and  using  the  primers 
(phoX  and  rotA)  from  Orchard  et  al.  (2009).  These  primers  were  designed  to  conserved 
regions  of  these  genes  and  have  been  used  to  amplify  a  broad  spectrum  of  Trichodesmium 
species.  Further,  the  expression  of  phoX  has  been  examined  in  a  representative  species 
from  a  second  Trichodesmium  clade  {Trichodesmium  spp.  H9-4)  (Chappell  and  Webb, 
2009)  and  was  found  to  have  the  same  expression  pattern  for  phoX  (data  not  shown). 
Melt  curves  for  each  reaction  were  performed  to  test  the  quality  of  the  product.  For  each 
sample  and  primer  pair  no  reverse  transcriptase  controls  and  no  template  controls  were 
run  to  test  for  DNA  contamination  and  specificity  of  amplification.  These  negative 
controls  did  not  typical  amplify,  but  in  the  rare  cases  where  there  was  any  amplification 
the  measured  template  concentration  was  at  least  two  orders  of  magnitude  less  than  the 
sample  template  concentration.  For  field  samples,  multiple  template  dilutions  were  tested 
for  each  sample  to  confirm  that  the  primers  had  the  same  efficiencies  as  they  did  in 
culture  (Orchard  et  ah,  2009),  and  therefore  confirm  that  the  reactions  were  not  inhibited. 
It  was  empirically  determined  that  when  the  template  was  diluted  at  least  1:5  for  all  the 
field  samples  tested  here,  the  reaction  was  not  inhibited. 

The  comparative,  or  AACt,  method  was  used  to  normalize  the  expression  of  each 
gene.  The  amplification  efficiency  and  validation  of  the  AACt  method  were  tested  for 
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each  primer  set  on  cDNA.  Samples  were  endogenously  normalized  to  the  constitutively 
expressed  rotA  (encoding  a  peptidyl-prolyl  cis-trans  isomerase)  (Orchard  et  ah,  2009). 
Fold  changes  were  calculated  using  the  Relative  Expression  Software  Tool  (REST), 
which  accounts  for  differences  in  efficiency  in  its  fold  change  calculation 
(http://www.gene-quantification.de/download.html).  In  both  field  and  culture  expression 
data  a  replicate  of  the  0.25  d'^  diluted  culture  (the  most  P  replete)  was  used  as  a  calibrator 
for  calculating  fold  changes.  The  p-values  were  calculated  using  a  pair-wise  fixed 
reallocation  randomization  test  (Pfaffl  et  ah,  2002). 

Results 

Semi-continuous  experiment 

Three  different  metrics  were  used  to  examine  P  physiology  in  the  semi- 
continuous  culture  experiment:  APA,  C:P  ratios  (using  both  total  cellular  P,  Ptc;  and 
intracellular  P,  Pic),  and  the  expression  of  phoX.  APA  normalized  to  chi  a  tracked 
inversely  with  dilution  rate  (Eigure  2C).  APA  from  cultures  growing  at  0.1  d'^  was 
significantly  higher  than  cultures  growing  at  0.25  d'^  (p<0.001).  The  ratio  of  C:Ptc 
ranged  from  96-156,  and  also  tracked  inversely  with  dilution  rate  (Eigure  2E).  This  range 
spans  from  slightly  below  the  Redfield  ratio  (106:1)  to  well  above.  Although  the  O.ld'^ 
and  0.15  d'^  diluted  cultures  had  similar  C:Ptc  ratios  of  149.7-I-/-0.524  and  156.1-1-/- 12. 45, 
the  C:Pic  for  the  0.10  d'^  cultures  was  much  higher  than  the  0.15  d'^  cultures.  The 
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cultures  diluted  at  0.1  d"'  had  a  significantly  higher  C:Ptc  and  C:Pic  than  the  cultures 
diluted  at  0.25  d'^  (p=0.04  and  p=0.004  respeetively). 

Like  C:P  ratio  and  APA,  phoX  expression  tracked  inversely  with  dilution  rate. 
The  expression  of  phoX  was  significantly  higher  (56  fold;  p>0.001)  in  the  0.1  d'*  diluted 
culture  compared  to  the  0.25  d"'  diluted  culture  (Figure  2D).  The  expression  of  phoX  in 
the  0.2  d'^  diluted  cultures  was  intermediate  with  a  fold  ehange  of  17  higher  than  the  0.25 
d'^  diluted  cultures.  Therefore,  the  expression  of  phoX  has  a  similar  relationship  as  the 
gene  it  eneodes  and  tracks  with  traditional  metrics  of  P  physiology. 

N2  fixation  rates  were  measured  for  each  culture.  The  cultures  diluted  at  0.25  d'^ 
had  the  highest  N2  fixation  rates  per  unit  ehl  a,  with  N2  fixation  decreasing  with 
increasing  P  limitation.  The  0.1  d"'  dilution  rate  cultures  had  a  significantly  lower  N2 
fixation  rate  than  the  eultures  diluted  at  0.25  d'*  (p<0.05)  (Figure  2A).  The  highest 
expression  for  nifH  in  the  0.25  d"'  diluted  sample  at  1.2  the  replicate  0.25  d'^  culture 
(Figure  2B).  The  nifH  gene  was  down  regulated  in  the  0.2  d'*  and  the  0.1  d'^  diluted 
culture  by  2.3  and  5.8  fold  respectively.  The  expression  of  nifH  was  significantly  lower 
(p<0.001)  for  the  0.1  d'^  diluted  eulture  relative  to  the  0.25  d'*  diluted  eulture.  Therefore, 
nifH  expression  showed  the  same  relative  pattern  between  the  cultures  as  N2  fixation. 

By  plotting  the  expression  of  phoX  versus  the  percent  reduction  in  N2  fixation 
(relative  to  the  maximal  N2  fixation  measured)  phoX  expression  can  be  calibrated  to  P 
driven  reductions  in  N2  fixation  (Figure  2F).  This  ealibration  ean  then  be  used  to  prediet 
reductions  in  N2  fixation  in  the  field  based  on  the  degree  of  phoX  expression.  The  data 
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was  divided  into  three  bins  indicating  the  level  of  phoX  expression  where  there  was  a  0- 
40%,  40-70%  or  above  70%  P  driven  reduction  in  N2  fixation. 

Basin-wide  comparisons  of  DIP,  APA  and  phoX  expression. 

Inorganic  phosphate  (DIP)  concentration  varied  with  region.  The  South  Pacific 
had  the  highest  DIP  concentration  ranging  from  43-114  nM,  and  the  North  Atlantic 
samples  had  significantly  lower  DIP  concentrations  ranging  from  0.5-17.6  nM  (Table  1; 
Figure  3)  (p<0.001).  Trichodesmium  APA  was  significantly  higher  in  the  North  Atlantic 
(0.23-0.56  nmol  P  h'^  colony'^)  than  the  South  Pacific,  (0.009-0.03  nmol  P  h'^  colony'*) 
(Table  1;  Figure  3)(p<0.001).  Trichodesmium  APA  also  trended  with  DIP  concentration, 
where  the  highest  APA  values  occurred  in  the  locations  with  the  lowest  DIP 
concentrations.  These  data  had  a  log-linear  relationship  with  (r  =0.83;  p>0.001).  While 
the  data  was  significant  correlated,  most  of  the  values  clustered  together  in  the  North 
Atlantic  with  low  DIP  and  high  APA,  or  in  the  South  Pacific  with  low  APA  and  higher 
DIP  (Figure  3C). 

The  expression  of  phoX  in  the  field  was  variable  depending  on  location.  The 
South  Pacific  samples  expression  ranged  from  5-17.5  fold  higher  than  the  0.25  d'*  diluted 
culture.  The  samples  from  the  North  Atlantic  were  higher  with  expression  ranging  from 
51-453  fold  higher  than  the  0.25  d'*  diluted  culture.  The  data  from  X0804  station  8  is 
included  in  the  graphs,  but  not  in  the  regression,  because  the  DIP  concentration  at  this 
station  is  anomalously  high.  The  expression  of  phoX  was  significantly  correlated  to  DIP 
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concentration  with  a  log-log  relationship  that  had  a  (r^=  0.75;  p>0.001)  (Figure  3 A  and 

B). 

Discussion 

Understanding  what  constrains  N2  fixation  is  critical  to  understanding  ocean 
biogeochemistry  and  the  limits  on  primary  productivity.  There  have  been  a  number  of 
studies  in  the  North  Atlantic  that  have  suggested  that  P  could  be  limiting  to  N2  fixation 
(Sanudo-Wilhelmy  et  ah,  2001;  Mills  et  ah,  2004;  Webb  et  ah,  2007)  and  to 
Trichodesmium  (Dyhrman  et  ah,  2002;  Mulholland  et  ah,  2002).  But  without  a 
quantitative,  genus-specific  metric  for  P  physiology,  it  is  difficult  to  assess  the  extent  of  P 
limitation  in  natural  populations  of  Trichodesmium  and  the  extent  to  which  P  could  be 
limiting  N2  fixation.  By  examining  the  expression  of  P-regulated  genes  in 
Trichodesmium  and  linking  those  expression  patterns  to  known  affects  on  N2  fixation 
rates  it  was  possible  to  identify  the  level  of  P  limitation  in  situ,  and  predict  reductions  in 
N2  fixation  due  to  P  limitation  in  Trichodesmium. 

Culture  experiments 

Trichodesmium  erythraeum  IMS  101  was  entrained  at  different  P-limited  growth 
rates  (by  varying  the  dilution  rate)  to  calibrate  the  expression  of  phoX  to  decreases  in  N2 
fixation  rate  due  to  P  limitation.  APA,  C:P  and  phoX  expression  tracked  inversely  with 
growth  rate,  N2  fixation,  and  nifH  expression.  The  cultures  with  the  highest  growth  rate, 
and  the  largest  input  of  new  P  per  day  had  the  lowest  APA  and  the  highest  N2  fixation 
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rate  when  normalized  to  chi  a.  Therefore,  N2  fixation  in  Trichodesmium  was  constrained 


by  P  supply.  In  fact,  the  N2  fixation  rate  in  the  cultures  with  the  lowest  growth  rate  (0.1 
d'^)  was  reduced  by  70%  compared  to  that  of  the  cultures  diluted  at  the  highest  growth 
rate  (0.25  d'^),  and  all  metrics  of  P  physiology  indicated  that  Trichodesmium  was  P 
limited  in  those  cultures.  A  similar  P  driven  reduction  in  N2  fixation  rate  of  66%  in 
Trichodesmium  semi-continuous  culture  was  observed  in  Hutchins  et  al  (2007).  This 
experiment  allows  for  the  calibration  of  phoX  expression  to  reductions  in  N2  fixation  due 
to  P  limitation  (Figure  2F),  which  can  then  be  used  to  predict  the  degree  to  which  N2 
fixation  is  limited  by  P  in  field  samples  by  measuring  phoX  expression. 

It  is  important  to  note  that  these  cultures  were  not  axenic,  and  both  C:P  ratios  and 
APA  were  influenced  by  heterotrophs  in  the  culture.  Additionally,  Trichodesmium 
colonies  in  the  field  are  rich  microbial  consortia,  and  APA  and  C:P  would  likewise  be 
affected  by  heterotrophs  associated  with  Trichodesmium  in  natural  populations. 
Furthermore,  APA  in  Trichodesmium  does  not  appear  to  be  actively  degraded  (Orchard  et 
ah,  2009),  and  therefore  the  APA  more  accurately  reflects  the  nutritional  history  of  the 
population  not  the  current  physiology.  So  while  APA  does  track  with  Trichodesmium 
physiology  this  not  the  ideal  metric  to  use  in  field  populations  to  identify  genus-specific  P 
limitation. 

The  C:P  ratio  is  frequently  used  as  a  measure  of  P  physiology  in  both  culture  and 
field  populations  (Falkowski,  1997;  Tyrrell,  1999).  However,  there  is  considerable 
plasticity  in  this  ratio  in  cyanobacteria  in  general,  and  in  Trichodesmium  in  particular, 
which  confounds  the  interpretation  (Bertilsson  et  ah,  2003;  Krauk  et  ah,  2006;  White  et 
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al.,  2006).  A  further  complication  is  that  a  portion  of  the  P  associated  with  the  cell 
appears  to  be  adsorbed  to  the  surface  (Sanudo-Wilhelmy  et  al.,  2004).  C:Pic  ratios 
followed  the  same  general  pattern  as  APA  and  the  inverse  of  N2  fixation  rates.  C:Ptc 
ratios  were  very  similar  for  cultures  grown  at  either  0.1  d'*  or  0.15  d'^  This  was 
unexpected  considering  the  APA  data,  which  suggests  that  the  cultures  diluted  at  0.10  d'^ 
were  much  more  limited  than  those  diluted  at  0.15  d'\  However,  based  on  C:Pic 
concentration  cultures  diluted  at  0.1  d'^  had  much  higher  C:P  ratios.  It  may  also  be  that 
oxalate  washing  removes  more  P  than  is  simply  adsorbed  to  the  surface,  or  that  with 
increased  APA  proportionally  more  P  was  adsorbed  to  the  outside  of  these  cells  as  a 
natural  step  before  being  assimilated  into  the  cell.  Taken  together,  C:Pic  tracks  best  with 
other  metrics  of  P  physiology,  but  given  the  variability  of  this  metric  observed  in  this 
study  and  elsewhere,  C:P  ratios  may  not  be  the  most  practical  method  for  inferring  P 
limitation  of  N2  fixation  rates  in  the  field. 

These  culture  results  demonstrate  that  the  expression  of  phoX  can  be  used  to 
identify  P  physiology  of  Trichodesmium  and  infer  affects  on  N2  fixation  by  this  genus. 
Utilizing  the  expression  of  phoX  has  the  added  benefit  of  being  genus-specific,  so  unlike 
APA  or  C:P,  examining  phoX  expression  will  not  be  affected  by  epibionts.  Therefore, 
using  phoX  expression  is  a  useful  metric  for  identifying  Trichodesmium  physiology  and 
identifying  areas  where  P  limits  N2  fixation. 

Basin-wide  comparisons  of  DIP,  APA  and  phoX  expression 
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The  expression  of  phoX  was  detectable  in  all  field  populations  of  Trichodesmium 
in  the  North  Atlantic  and  South  Pacific.  To  our  knowledge  nifH  is  the  only  other 
Trichodesmium  functional  gene  whose  expression  has  been  measured  quantitatively  in  a 
marine  system  (Church  et  ah,  2005;  Foster  et  ah,  2009).  The  expression  of  phoX  in  the 
field  shows  a  log-log  relationship  with  DIP  concentration.  It  is  important  to  note  that 
phoX  expression  is  regulated  by  P  physiology  not  DIP  supply.  Therefore,  whether 
Trichodesmium  is  utilizing  DIP,  or  DOP,  as  a  P  source,  phoX  is  down-regulated.  This 
has  been  demonstrated  in  Synechococcus  where  phoX  was  not  highly  expressed  when 
growing  on  phosphonate  (a  DOP  compound)  as  compared  to  cultures  that  were  P-starved 
(Adams  et  ah,  2008).  Additionally  spiking  either  P  ester  or  phosphonate  to  a  P-starved 
culture  of  Trichodesmium  results  in  the  down-regulation  of  phoX  within  4  hours  (data  not 
shown).  Therefore,  this  tight  correlation  between  phoX  and  DIP  implies  that  the  bio- 
available  P  pool  is  probably  correlated  to  the  DIP  concentration. 

There  was  also  a  significant  relationship  between  APA  and  DIP.  Although  the 
majority  of  the  samples  clustered  together  at  very  high  rates  for  the  North  Atlantic 
samples  and  low  rates  for  the  South  Pacific  samples.  Given  that  APA  is  not  a 
Trichodesmium  specific  signal  and  that  previous  culture  results  suggest  that  APA  may  not 
be  actively  degraded,  we  would  expect  to  see  generally  higher  values  in  the  North 
Atlantic. 

Within  the  North  Atlantic  there  were  detectable  differences  in  the  expression  of 
phoX  and  therefore  gradations  in  the  degree  of  P  limitation  in  Trichodesmium  within  this 
ocean  gyre.  Further  demonstrating  that  phoX  is  a  very  sensitive  metric  of  P  physiology. 
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It  is  striking  that  even  in  the  very  low  P  North  Atlantic  small  changes  in  DIP 
concentration  can  led  to  large  changes  in  phoX  expression,  suggesting  even  more 
variation  in  the  degree  of  P  limitation  in  this  area  than  previously  considered. 

The  high  level  of  expression  of  phoX  in  the  North  Atlantic  Trichodesmium 
samples  implies  that  Trichodesmium  growth  rate  and  N2  fixation  by  were  depressed,  and 
C:P  ratios  was  elevated.  This  is  consistent  with  previous  studies  from  the  North  Atlantic 
that  have  shown  elevated  C:P  ratios  for  Trichodesmium  (Krauk  et  ah,  2006;  White  et  ah, 
2006).  Based  on  the  level  of  phoX  expression,  Trichodesmium  N2  fixation  is  P-limited  in 
the  North  Atlantic  by  over  40%  and  in  many  cases  over  70%  relative  to  maximal 
potential  rates  in  culture.  These  data  suggest  that  if  Trichodesmium  were  re-fed  P  in  the 
North  Atlantic,  then  N2  fixation  rate  would  increase  in  this  region  by  as  much  as  70% 
(assuming  no  co-limitation  is  occurring). 

P  limitation  in  other  gyres 

There  was  little  change  in  phoX  expression  above  7  nM,  but  below  7  nM  phoX 
expression  increased  rapidly  (Figure  3A).  This  implies  that  there  could  be  a  threshold 
concentration  of  DIP  (representing  bio-available  P)  where  P  becomes  limiting.  Moutin  et 
al  (2005)  calculated  a  similar  threshold  DIP  concentration  of  9  nM  for  Trichodesmium  in 
the  South  Pacific  based  on  DIP  uptake  kinetics.  Assuming  that  this  threshold  DIP 
concentration  is  applicable  for  other  regions  of  the  ocean  (in  addition  to  the  South  Pacific 
and  North  Atlantic),  it  can  be  used  to  predict  where  P  is  limiting  to  Trichodesmium  in  the 
ocean.  The  western  North  Atlantic  has  chronically  low  DIP  concentration  in  the 
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subnamolar  to  10  nM  range  (Wu  et  al.,  2000;  Lomas  et  al.,  2009),  therefore  we  would 
predict  that  P  general  limits  Trichodesmium  growth  and  N2  fixation  in  this  region.  DIP 
concentrations  in  the  North  Pacific  are  much  more  variable  ranging  from  6-100  nM 
(Bjorkman  et  al.,  2000;  Grabowski  et  al.,  2008),  therefore  Trichodesmium  could  at  times 
be  experiencing  P  limitation  in  the  North  Pacific  but  this  would  be  transient  compared  to 
the  North  Atlantic.  While  DIP  concentrations  in  the  South  Atlantic  and  the  South  Pacific 
tend  to  be  higher  than  the  7  nM  threshold,  there  are  locations,  at  least  for  the  South 
Pacific  where  the  DIP  concentration  has  dropped  below  this  threshold  (Moutin  et  al., 
2005).  This  suggests  that  while  P  is  probably  not  limiting  for  most  of  the  year,  there  are 
times  where  DIP  could  reach  limiting  concentrations.  In  the  Mediterranean  P 
concentrations  can  be  sub-nanomolar  (Thingstad  et  al.,  2005),  and  if  the  Trichodesmium 
in  this  region  have  a  similar  P  threshold  then  P  would  be  limiting.  As  this  study  would 
predict  based  on  threshold  P  concentrations,  during  a  large-scale  P  addition  experiment 
that  was  conducted  in  the  Mediterranean,  there  was  an  increase  in  N2  fixation  rate  relative 
to  outside  of  the  P  addition  area  (Rees  et  al.,  2006). 

It  has  been  suggested  that  diazotrophy  is  enhanced  under  higher  CO2  conditions 
(Hutchins  et  al.,  2007).  Increases  in  N2  fixation  could  drive  regimes  towards  P  limitation 
(Hutchins  et  al.,  2007).  It  has  further  been  suggested  that  increased  N2  fixation  is  driving 
the  North  Pacific,  a  traditionally  N-limited  system,  towards  P  limitation  (Karl,  1999),  and 
increased  N2  fixation  would  only  further  drive  the  system  towards  P  limitation. 
Therefore,  while  currently  there  are  few  places  in  the  ocean  where  N2  fixation  is 
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predicted  to  be  limited  by  P  (Moore  et  al.,  2004),  P  limitation  could  become  much  more 
wide  spread  in  the  future. 

Conclusion 

The  expression  of  phoX  is  a  useful  quantitative  species-specific  metric  for  P 
physiology  in  Trichodesmium.  Using  a  semi-continuous  culture  approach  we  determined 
levels  of  phoX  indicative  decreased  N2  fixation  due  to  P  limitation.  In  field  populations 
there  is  a  tight  relationship  between  DIP  concentration  and  phoX  expression.  Further, 
these  data  suggest  that  the  North  Atlantic  is  chronically  P-limited,  but  that  there  are 
detectable  differences  in  the  degree  of  limitation  using  phoX  expression.  All  the  samples 
measured  in  the  North  Atlantic  had  a  level  of  phoX  expression  greater  or  very  near  the 
level  at  which  there  was  a  70%  reduction  in  N2  fixation  rate  in  culture,  but  South  Pacific 
samples  were  not  P-limited  based  on  these  metrics.  The  relationship  between  phoX 
expression  and  DIP  concentration  suggests  that  there  is  a  threshold  concentration  of  bio- 
available  P  that  is  proportional  to  7  nM  DIP,  and  at  this  threshold  P  starts  to  become 
limiting  for  N2  fixation  in  Trichodesmium.  While  there  are  few  areas  of  the  ocean  that 
have  DIP  concentrations  consistently  below  this  threshold,  seasonal  variability  can  lead 
to  periods  of  P  limitation.  Additionally,  if  CO2  concentrations  continue  to  rise,  and 
oligotrophic  systems  are  driven  towards  lower  DIP  concentration,  then  P  limitation  of  N2 
fixation  may  become  more  common.  This  work  utilized  molecular  indicators  of 
physiology  and  link  molecular  biology  to  biogeochemical  processes.  These  data 
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underscore  the  potential  for  combining  these  techniques  to  probe  genus  or  species- 
specific  responses  to  environmental  stimuli. 
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Table  1.  Station  locations,  sampling  times,  DIP  concentrations,  Trichodesmium  consortia  alkaline  phosphatase  activity,  and 
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Figure  1.  Station  locations  from  (A)  the  South  Pacific  and  (B)  the  North  Atlantic.  In  (B) 
open  circles  indicate  stations  from  May  2008  (X0804),  the  black  circle  indicates  a  station 
from  June  2007  (X0705)  and  the  grey  circle  indicates  a  station  from  October  2007 
(BVAL39). 
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Figure  2.  Trichodesmium  N2  fixation  rate  (A),  nifH  expression  (B),  alkaline  phosphatase  activity 
(APA)  (C),  phoX  expression  (D),  C:P  ratio  (E)  and  the  relationship  of  phoX  expression  and  % 
reduction  in  N2  fixation  relative  to  maximal  measured  rates  (F)  from  semi-continuous  culture 
experiments.  In  plot  F  data  is  hinned  to  show  the  predicted  reduction  in  N2  fixation  rate  at 
different  levels  of  phoX  expression.  Plots  A,  C  and  E  represent  the  average  value  of  duplicate 
bottles,  and  error  bars  indicate  standard  deviations.  For  plots  B,  D  and  F  error  bars  indicate  the 
standard  error  of  the  fold  change  of  triplicate  measurements  on  a  single  biological  replicate 
between  the  sample  expression  and  the  calibrator  expression  (a  replicate  of  the  0.25  d'^  diluted 
culture).  The  (*)  indicates  a  p-value  of  less  than  0.05  for  the  sample  relative  to  the  0.25  d  ' 
diluted  culture  using  a  t-test  (A,  C  and  E),  or  a  pair-wise  fixed  reallocation  randomization 
analysis  (B  and  D).  In  plot  E,  molar  C:P  ratios  from  cultures  are  shown  for  both  total  cellular 
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Figure  3.  Plots  of  phoX  expression  versus  DIP  concentration  on  a  linear  (A)  and  log-log 
scale  (B),  alkaline  phosphatase  activity  (APA)  versus  DIP  concentration  (C).  The  dashed 
line  in  A  indicates  the  7  nM  DIP  threshold  line.  The  dashed  line  in  B  indicates  the  fold 
change  for  phoX  in  the  semi-continuous  experiment  where  N2  fixation  was  reduced  by 
70%  relative  to  the  maximum  for  the  semi-continuous  cultures.  Error  bars  in  plots  A  and 
B  represent  the  standard  error  of  the  fold  change.  The  solid  lines  indicate  the  regression 
lines.  A  power  regression  for  plots  A  and  B  (r  =0.75)  and  an  exponential  regression  was 
used  to  fit  plot  C  (r^=0.83). 
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Chapter  Six 


Summary  and  Future  Directions 


Summary 

The  marine  cyanobacterium  Trichodesmium  is  a  major  contributor  to  dinitrogen 
(N2)  fixation  and  primary  production  in  tropical  and  subtropical  oceans.  These  regions 
are  often  characterized  by  low  phosphorus  (P)  concentrations,  and  P  starvation  may  limit 
growth  and  N2  fixation  in  Trichdoesmium.  This  thesis  focused  on  examining  three  major 
components  of  P  metabolism:  dissolved  inorganic  P  (P,j  uptake  (chapter  two  and  four), 
dissolved  organic  P  (DOP)  hydrolysis  (chapter  two  and  four)  and  polyphosphate  (polyP) 
formation  and  degradation  (chapter  three).  Additionally,  this  thesis  examines  the  extent 
of  P  limitation  and  its  possible  affect  on  N2  fixation  in  Trichodesmium  from  the  North 
Atlantic  and  South  Pacific  (chapter  five). 

Recent  studies  have  hypothesized  that  between  12-60%  of  the  primary 
productivity  in  the  Sargasso  Sea  is  supported  by  P  ester  not  Fi  (Mather  et  ah,  2008; 
Lomas  et  ah,  2009).  Chapter  two  addresses  the  relative  importance  of  P,  and  P  ester  in 
supporting  the  growth  of  Trichodesmium,  the  microbial  community,  and  other 
taxonomically  significant  groups  from  the  literature  (Casey  et  ah,  2009).  The  in  situ 
uptake  rates  suggest  that  the  contribution  of  P  ester  to  total  P  uptake  can  be  greater  than 
25%  for  Trichodesmium,  the  microbial  community,  and  key  phytoplankton  groups, 
depending  on  the  station.  The  P  ester  contribution  to  Trichodesmium  is  over  25%  at  half 
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of  the  stations  examined.  Further  work  in  other  systems  with  variable  P  chemistry  will 
indicate  whether  this  frequency  is  similar  in  higher  P  systems,  or  to  what  extent  the  P 
ester  uptake  rate  changes  with  ambient  DOP  over  a  larger  dataset.  A  comparison  of  the 

OO  TO 

Km  values  for  both  P;  and  P  ester  uptake  suggests  that  Trichodesmium  is  a  better 
competitor  for  P  ester  than  for  P,  relative  to  the  microbial  community.  This  is  consistent 
with  previous  work  (Sohm  and  Capone,  2006),  which  focused  solely  on  P,  uptake 
kinetics.  Additionally,  this  chapter  focused  on  enzyme  activity  and  uptake  assays,  and  it 
highlights  the  importance  of  pairing  tracer  level  experiments  (which  measure  in  situ 
uptake  rates)  to  experiments  where  saturating  concentrations  of  substrate  are  added 
(which  measure  the  Vmax  of  enzyme  hydrolysis).  Trichodesmium  growth  rates  calculated 

TO  oo 

from  the  P,  and  P  ester  uptake  parameters  indicate  that  only  a  small  fraction  of  the  P 
ester  pool  needs  to  be  bio-available  to  meet  P  demand,  which  further  supports  the 
importance  of  P  ester  to  this  group  in  the  Sargasso  Sea.  Taken  together,  these  data 
underscore  the  importance  of  P  esters  to  supporting  production  by  Trichodesmium, 
picocyanobacteria,  and  small  eukaryotes  in  the  Sargasso  Sea. 

P  in  the  form  of  polyP  has  recently  been  detected  in  coastal  systems  (Diaz  et  ah, 
2008)  and  has  previously  been  observed  in  Trichodesmium  (Romans  et  ah,  1994),  but 
polyP  has  never  been  examined  in  a  very  low  P  environment  such  as  the  Sargasso  Sea. 
Chapter  three  shows  that  polyP  can  be  a  major  component  of  Trichodesmium  cellular  P  in 
the  Sargasso  Sea,  despite  other  metrics  of  P  physiology  indicating  that  these  populations 
were  P-limited.  While  it  is  unclear  what  exact  factors  induce  the  formation  of  polyP  in 
this  system,  the  data  are  consistent  with  a  possible  overplus-type  response.  The  over-plus 
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response  occurs  when  P-limited  cells  accumulate  polyP  after  being  refed  P  as  opposed  to 
luxury  uptake  where  P  replete  cells  take  up  P  above  their  nutritional  demand  for  P  and 
accumulate  that  P  as  polyP.  It  may  be  that  Trichodesmium  is  taking  advantage  of  fluxes 
in  P  supply,  and  that  polyP  is  a  physiological  fingerprint  of  this  variability  in  the  Sargasso 
Sea.  Additionally,  if  polyP  formation  is  a  common  trait  in  phytoplankton,  these  data 
suggest  that  polyP  might  be  an  important  additional  bio-available  component  of  DOP. 
This  study  highlights  the  potential  importance  of  polyP,  and  suggests  the  further  studies 
into  the  dynamics  of  polyP  formation  in  oligotrophic  systems  are  a  critical  area  of  future 
research. 

Chapter  four  studies  the  affect  of  P  physiology  on  alkaline  phosphatase  activity 
(APA)  and  gene  expression  patterns  in  Trichodesmium.  Two  copies  of  the  high  affinity 
phosphate  binding  protein  {pstS  and  sphX)  and  two  putative  alkaline  phosphatase  (phoA, 
and  phoX)  genes  are  present  in  all  examined  species  of  Trichodesmium.  Additionally,  the 
expression  of  phoA,  phoX,  and  sphX  are  all  responsive  to  variations  in  P  supply,  while 
pstS  is  not.  These  data  suggest  that  the  expression  of  sphX,  phoA  and  phoX  could  provide 
a  useful  marker  for  P  limitation  in  natural  populations.  Combining  these  measurements 
with  APA  would  also  provide  insights  into  the  nutritional  history  of  this  organism.  This 
is  the  first  demonstration  of  phoX  P-regulation  in  a  marine  cyanobacterium.  The 
presence  of  a  P-regulated  phoX  in  Trichodesmium  could  provide  a  means  to  access  the  P 
ester  pool  without  an  increased  zinc  (Zn)  requirement  (Wu  et  ah,  2007).  In  oligotrophic 
regimes,  where  P  esters  are  the  dominant  P  pool  and  Zn  concentrations  are  low  (Jakuba  et 
ah,  2008),  having  a  PhoX  type  alkaline  phosphatase  could  allow  Trichodesmium  to 
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occupy  a  unique  niche  relative  to  cyanobacteria  solely  using  a  PhoA  type  alkaline 
phosphatase,  which  typically  requires  Zn  (Shaked  et  ah,  2006).  Molecular 
characterization  of  responses  to  P  limitation  lay  a  foundation  for  better  understanding  the 
dynamics  of  Trichodesmium  P  limitation  in  the  sea. 

Chapter  five  calibrates  the  expression  of  phoX  to  decreases  in  growth  and  N2 
fixation  due  to  P  limitation  in  Trichodesmium.  The  expression  of  phoX  is  then  used  to 
examine  the  extent  of  P  limitation  in  situ.  Trichodesmium  was  P-limited  in  the  North 
Atlantic,  and  the  expression  of  phoX  is  indicative  of  as  much  as  a  70%  reduction  in  N2 
fixation  by  this  genus,  but  P  was  no  limiting  Trichodesmium  in  the  South  Pacific  from  the 
transect  examined  here.  These  data  also  suggest  that  P  becomes  limiting  to 
Trichodesmium  at  approximately  7  nM  Pi,  which  agrees  with  previous  estimates  of  a 
threshold  concentration  for  Trichodesmium  of  9  nM  (Moutin  et  ah,  2005).  At  this 
concentration  P  would  not  be  limiting  to  Trichodesmium  in  the  South  Atlantic,  but  could 
be  limiting,  at  times,  in  the  North  Pacific.  These  results  support  models  predicting  P 
limitation  of  N2  fixation  the  North  Atlantic  but  not  South  Pacific  or  South  Atlantic 
(Moore  et  ah,  2004).  While  N2  fixation  by  Trichodesmium  is  higher  in  the  North  Atlantic 
than  many  other  oligotrophic  gyres,  probably  due  to  iron  (Fe)  stress  (Kustka  et  ah,  2002; 
Moore  et  ah,  2004)  in  other  regions,  our  results  suggest  that  N2  fixation  by  this  genus 
could  be  significantly  higher  with  an  increased  P  supply  to  the  North  Atlantic. 

Future  directions 
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As  is  typical,  this  thesis  introduces  just  as  many  questions  and  future  avenues  of 
research  as  it  addresses.  For  example,  chapter  two  highlights  the  importance  of  P  ester  to 
support  primary  production  in  Trichodesmium,  the  microbial  community,  and  a  few  other 
taxonomically  important  genera.  Future  work  examining  the  importance  of  phosphonate 
will  allow  us  to  understand  the  relative  importance  of  different  DOP  bond  classes. 
Extending  this  work  to  include  more  taxonomic  groups  and  more  ocean  gyres  with 
differing  P  chemistry  will  help  to  determine  the  global  significance  of  DOP  as  a  P  source. 

The  dynamics  of  polyP  formation  and  degradation  remains  somewhat  enigmatic. 
While  polyP  is  clearly  important  to  Trichodesmium  its  prevalence  in  other  oligotrophic 
microbes  is  unknown.  This  compound  has  the  potential  to  be  very  important  to  the  P 
cycle  in  the  ocean  gyres,  and  its  synthesis  and  degradation  rates  should  be  assayed  and 
compared  to  standing  stock  concentrations  to  better  constrain  the  flux  of  polyP.  Future 
culture  work  studying  the  conditions  that  promote  polyP  formation  are  also  needed  to 
determine  the  role  this  compound  plays  in  cells. 

Chapter  four  examines  key  components  of  the  Pho  regulon  in  Trichodesmium,  but 
there  are  many  more  genes  to  examine.  For  example,  Trichodesmium  has  genes  encoding 
enzymes  involved  in  phosphonate  scavenging  whose  regulation  have  only  been  examined 
qualitatively  (Dyhrman  et  ah,  2006).  Future  work  characterizing  the  regulation 
phosphonate  scavenging  genes  will  broaden  our  understanding  of  DOP  hydrolysis  in 
Trichodesmium.  Furthermore,  the  regulation  of  the  genes  examined  in  this  study  suggest 
that  they  are  correctly  annotated,  but  to  confirm  these  annotations  the  proteins  need  to  be 
characterized.  This  is  particularly  important  for  the  two  alkaline  phosphatases:  PhoA  and 
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PhoX.  Determining  the  substrate  specificity,  cellular  localization  and  metal  requirements 
of  PhoX  and  PhoA,  will  allow  for  us  to  understand  which  P  esters  are  bio-available,  and 
constrain  trace  metal  requirements.  Considering  that  Zn  may  also  often  be  limiting  to 
growth  (Shaked  et  ah,  2006;  Jakuba  et  ah,  2008),  the  presence  of  the  PhoX  enzyme 
suggests  an  interesting  adaptation  to  avoid  Zn-P  co-limitation.  Therefore,  additional 
studies  to  determine  growth  rates  and  gene  expression  patterns  over  ranges  in  Zn  and  P 
concentrations  are  warranted. 

Utilizing  gene  expression  as  a  marker  of  P  limitation  in  Trichodesmium  allowed 
for  a  unique  and  powerful  view  into  the  physiological  ecology  of  this  biogeochemically 
significant  organism  in  the  North  Atlantic  and  South  Pacific.  The  screening  for  more 
Trichodesmium  populations  across  a  broader  range  of  Fi  concentrations  and  oligotrophic 
systems  will  further  elucidate  the  relationship  between  Pz  concentration  and 
Trichodesmium  physiology.  For  example,  will  the  7  nM  threshold  for  P  limitation  be  the 
same  for  Trichodesmium  from  the  South  Atlantic,  or  will  these  populations  have  a  higher 
Pz  threshold?  Also,  future  studies  can  utilize  additional  genes  to  examine  the  importance 
of  multiple  stressor  simultaneously.  Utilizing  the  expression  of  genes  that  are  diagnostic 
of  cellular  physiology  allows  for  the  simultaneous  analysis  of  multiple  metabolism 
pathways.  For  example,  idiA,  is  regulated  by  Fe  availability  (Chappell  and  Webb,  2009). 
Examining  the  regulation  of  both  phoX  and  idiA  from  the  same  population  would  allow  a 
clearer  understanding  of  the  relationship  between  P  and  Fe,  and  to  explore  the  possibility 
of  co-limitation.  Furthermore,  developing  this  type  of  assay  for  other  taxa  will  provide 
new  insights  into  what  limits  primary  production  in  the  ocean. 
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These  areas  of  research  will  expand  our  knowledge  of  Trichodesmium  P 


physiology,  the  importance  of  DOP  in  driving  primary  production,  and  the  role  nutrients 
play  in  controlling  N2  fixation. 
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